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Abstract. The Koch snowfiake KS is a nondifferentiable curve. Hence, any 
attempt to define reflection in the boundary may seem like an exercise in 
futility. In this paper, for each integer n > 0, we describe the periodic orbits of 
the prefractal billiard Q(KS n ) (the nth inner rational polygonal approximation 
of the Koch snowfiake billiard). Moreover, we use this information in order 
to define and describe a particular collection of periodic orbits of the Koch 
snowfiake billiard Q(KS). 

In the finite case, an orbit of fl(KS„) can be reduced to its Poincare 
section, which simply amounts to a finite collection of points in the boundary of 
the prefractal billiard. We show that, for each n > 0, the collection of directions 
for which the billiard flow on Q(KS n ) is closed is exactly the collection of 
directions for which the billiard flow on Q(KSo) is closed. Such a result relies 
on the fact that the corresponding flat surface S(KS n ) (n > 1) is shown to be 
a branched cover of the flat surface S(KSo), the hexagonal torus. Extending 
this result, we define what we call a compatible sequence of periodic orbits. 
Focusing on the direction given by an initial angle of rr/3, we define 1) a 
compatible sequence of piecewise Fagnano orbits, 2) an eventually constant 
compatible sequence of orbits and 3) a compatible sequence of generalized 
piecewise Fagnano orbits. 

In the case of the infinite (fractal) billiard table, we will describe what 
we call stabilizing periodic orbits of the Koch snowfiake billiard Q(KS). An 
eventually constant compatible sequence of periodic orbits is comprised (for all 
but finitely many) of ^-orbits (or what we also call stabilizing periodic orbits). 
We show that the trivial limit of an eventually constant compatible sequence 
of periodic orbits is, in fact, a periodic orbit of Q(KS). In a sense, we show 
that it is possible to define billiard dynamics on a Cantor set. 

In addition, we will discuss the geometric and topological properties of 
what we call the footprint of a piecewise Fagnano orbit. We will show that the 
inverse limit of the footprints of orbits of the prefractal approximations (or, 
the Poincare sections of the respective orbits) exists in a specific situation and 
provide a plausibility argument as to why such an inverse limit of footprints 
should constitute the footprint of a well-defined periodic orbit of Q(KS). Us- 
ing, in particular, known results for the inverse limit of a sequence of finite 
spaces, we deduce that the footprint (i.e., the intersection of the orbit with the 
boundary) of a piecewise Fagnano orbit is a topological Cantor set and even, 
a self-similar Cantor set. 

We allude to a possible characterization of orbits with an initial direction 
of 7r/3. That is, we provide support for a complete description of periodic orbits 
in the direction of it/ 3. Such a characterization would allow one to describe 
an orbit with an initial direction of n/3 of the Koch snowfiake billiard as 
either a piecewise Fagnano orbit, a stabilizing orbit or a generalized piecewise 
Fagnano orbit. We then close the paper by discussing several outstanding 
open problems and conjectures about the Koch snowfiake billiard Q(KS), the 
associated ' fractal flat surface' and possible connections with the associated 
fractal drum T>(KS) via fractal analogues of Gutzwiller-like trace formulae. 

These problems and conjectures have natural counterparts for other frac- 
tal billiards. In the long-term, the present work may help lay the foundations 
for a general theory of fractal billiards. 



1. Introduction 

The Koch snowfiake curve, as depicted in Figure HJ is a fractal. In particu- 
lar, it is the union of three self-similar Koch curves, with the Koch curve being a 
continuous, nowhere differentiable curve with infinite length (see Figures [2] and [3j) . 
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Consequently, any attempt to construct a line tangent to the Koch snowflake curve 
may seem like an exercise in futility. This poses a unique problem for defining the 
trajectory of a billiard ball (i.e., a pointmass traversing the interior of the planar 
region bounded by the Koch snowflake K S). Specifically, when this pointmass col- 
lides with the boundary K S with unit speed, the absence of a well-defined tangent 
results in multiple choices for the angle of reflection, meaning there is a priori no 
well-defined angle of reflection. In [LaNiel] . we provided experimental evidence 
in support of the existence of certain periodic orbits of the Koch snowflake billiard 
il(KS). Moreover, in [LaNiel] . we stated several conjectures about the existence 
of a well-defined billiard Q(KS) and the dynamical equivalence between the con- 
jectured billiard flow on tl(KS) and the associated geodesic flow on the proposed 
corresponding 1 fractal flat surface' . One of the main objectives of the present paper 
is to investigate what one means by reflection in the snowflake boundary and to 
establish the existence and describe the topological and geometric properties of par- 
ticular families of periodic orbits (and/or of their footprints) of the Koch snowflake 
billiard Vt{KS). 

In short, the point of view adopted here is to define certain "periodic orbits" 
of fl(KS) as suitable (inverse) limits of certain "compatible sequences" of periodic 
orbits of its (inner) rational polygonal billiard approximations Q(KS n ). Using this 
definition and a study (conducted in §3) of the periodic orbits of the nth prefractal 
billiard approximation Q(KS n ), for each fixed n > 0, we characterize and describe 
(in terms of the ternary expansion of their initial basepoint) the periodic orbits 
with an initial direction making an angle of 7r/3 with the horizontal in D,(KSq). 

More specifically, we are able to construct what we call the footprint of the 
primary piecewise Fagnano and piecewise Fagnano orbits, and stabilizing periodic 
orbits of il(KS) (all in the direction of it/3). As of now, the only family of well- 
defined orbits of the Koch snowflake billiard ft(KS) is the family consisting of what 
we call stabilizing orbits. That which we propose to be a piecewise Fagnano orbit 
of Sl(KS) has a footprint J-(xq) that is the inverse limit of footprints (of piece- 
wise Fagnano orbits) of the prefractal approximations. A footprint of a prefractal 
approximation amounts to the Poincare section of the billiard map describing the 
billiard flow on the corresponding phase space. While we say "piecewise Fagnano 
orbit," we are making an abuse of language in that we do not mean to imply that 
an orbit actually exists, but that whatever the orbit truly is, it has a footprint 
J-(xq). Furthermore, even less is known about what we have called the generalized 
piecewise Fagnano orbits. Again, there really is no orbit to speak of, nor is there 
any footprint to speak of. We discuss all of these 'orbits' in §4 and §5 with differing 
degrees of rigor, sometimes only providing a plausibility argument as to why a given 
orbit should have a particular property. 

This paper draws upon various subjects in mathematics. So as to accommo- 
date a diverse audience of readers, we make a considerable effort in developing the 
necessary background material in §2. In particular, we give a thorough description 
of the billiard flow associated with a billiard table £l(B) with (piecewise) smooth 
boundary B. We also recall the notion of a flat surface and how one can construct 
a flat surface from a rational polygonal billiard table (that is, a planar billiard table 
whose boundary is a polygon with interior angles that are rational multiples of 7r) . 
In this context, a flat surface is a mathematical device used to rigorously describe 
the billiard flow on il(B) in terms of the geodesic flow on the surface. In addition, 
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Figure 1. The Koch snowflake curve K S and its prefractal ap- 
proximations KS n , for n = 0, 1, 2, 6. 






Figure 2. The Koch curve KC and its prefractal approximations 
KC n , for n = 0, 1, 2, 6. The Koch curve is self-similar, meaning 
that there are scaled (and rotated/reflected) copies of the Koch 
curve found as subsets of KC. In this figure, we circle four copies 
of KC scaled by 1/3. In general, one can find 4™ copies of the 
Koch curve scaled by 1/3™ as subsets of KC such that the disjoint 
union (disjoint except at the endpoints of each copy) comprises the 
whole Koch curve. This is, in fact, the essence of self-similarity. 
By abuse of language, we say that the Koch snowflake curve itself 
is "self-similar" ; see Figure [31 



Figure 3. The self-similarity of the Koch snowflake curve KS. 
Shown here is the Koch snowflake curve KS, viewed as the union 
of three isometric, abutting copies of the Koch curve. 

we recall the definition of inverse limit and explain how the Cantor set can be 
viewed as the inverse limit of an inverse limit sequence of its prefractal approxima- 
tions, denoted by with the index n corresponding to the approximation with 
2™ many points. 

Similarly, the snowflake curve KS can be viewed as the inverse limit of its 
prefractal approximations KS n . Here, KS n is the nth (inner) polygonal approx- 
imation to KS, and hence defines a rational polygonal billiard table il(KS n ). 
Since the theory of rational polygonal billiards is very well developed (see, e.g., 
[GaStVo,Gul,GuJul-2,HuSc,KaHa2,KaZe,Mas,MasTa,Vel-3,Vo,Zo]), it is then nat 
ural to define the dynamics on the fractal "billiard table" Sl(KS) in terms of the 
dynamics on its prefractal approximations £l(KS n )- As a result, much of the focus 
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of this paper will be to first obtain a good understanding of the periodic orbits of 
fl(KSn), and then to provide a plausibility argument as to why we can view piece- 
wise Fagnano orbits on fl(KS) as suitable (inverse) limits of appropriate sequences 
of piecewise Fagnano orbits on Cl(KS n ), for n > 0. 

It is our intention that those familiar with the topic of mathematical billiards 
and not fractal geometry find readily accessible the basic notions of self-similarity 
and that of an iterated function system (IFS). We briefly describe these notions 
by means of a simple example in i j4.ll and refer the reader to various references 
for further details (e.g., |Bal IEd[ |Fa| ) . So as to accommodate readers from the 
physical sciences, we also attempt to explain the necessary concepts from topology 
and geometry. As such, the interested reader will find references to [Maj for further 
details on covering spaces, [McLj for category theory, [Bp] for general topology and 
[HoYoj for the specialized topic of inverse and direct limits in the context of the 
category of topological spaces. 

In §3, we begin the discussion of our results. We prove that the prefractal flat 
surface S(KS n ) associated with the Koch snowflake prefractal billiard Q(KS n ) is 
a branched cove$\ of the hexagonal torus S(KSq). We then use this result to show 
that initial directions of periodic orbits in the billiard Q(KSq) are exactly the initial 
directions of periodic orbits in the billiard Q(KS n ), and vice- versa. 

This key fact serves as the foundation for §4. It is there that we develop much 
of the machinery necessary to describe what we call compatible sequences of periodic 
orbits. Further discussion about the nature of particular points in the unit interval 
/ (which we always view as the base of the equilateral triangle KSo := A) gives rise 
to specific compatible sequences. We note that §4 is very dense, serving as a strong 
foundation for §5. While the addressing system used in §5 (and introduced and 
used in [LaPaj ) may seem to make some of the tools developed in §4 redundant, in 
principle, many of the proofs of the results in §5 demonstrate the interconnectedness 
of the two sections, making the implicit (and even explicit) dependence of §5 on §4 
readily apparent. 

Therefore, we model the structure of §5 on that of §4, so as to allow the reader to 
draw parallels more quickly and to see how the results in §4 influence our later study 
of the periodic orbits of the Koch snowflake billiard. In §4, we show that directions 
for which an orbit & n {x® n , 0°) of Q(KS n ) are periodic are exactly the same for which 
an orbit &o(xq, 9q) of Q(KSo) is periodic. This aids us in constructing what we call 
a compatible sequence of closed orbits. We focus our investigation on orbits with an 
initial direction of tt/3. As such, we describe what we call a compatible sequence of 
piecewise Fagnano orbits, an eventually constant compatible sequence of orbits and 
compatible sequence of generalized piecewise Fagnano orbits. The period and length 
of piecewise Fagnano orbits, c € -orbits and generalized piecewise Fagnano orbits of 
fl(KS n ) are given in terms of the ternary representation of the initial basepoint 
of the initial orbit of the respective compatible sequence of periodic orbits. In 
£15.3.11 and £15.41 we describe the topological and geometric properties of what we 
call the footprint of a piecewise Fagnano orbit and stabilizing orbits (or ^-orbits), 
respectively. In §5.5[ we then provide a plausibility argument for the existence 
of what we call a piecewise Fagnano orbit of the Koch snowflake billiard Q(KS). 
Finally, we close §5 by conjecturing the existence of generalized piecewise Fagnano 
orbits of the Koch snowflake billiard fl(KS). 



We briefly discuss the definition of branched cover in §3. 



6 



MICHEL L. LAPIDUS AND ROBERT G. NIEMEYER 



Considering the fact that the field of "fractal billiards" is still in its infancy, we 
provide many open questions and conjectures in §6. In particular, we stress that 
tt(KS) does not constitute a well-defined mathematical billiard, in the sense that 
we have not provided a well-defined phase space, let alone a geodesic flow on such 
a phase space. Such a mathematical object has yet to be precisely defined, but the 
work we have laid out in §5 and the remarks made in §6 indicate a possible path 
for constructing such a phase space and geodesic flow. 

In addition to determining the nature of such a geodesic flow and whether it 
could be dynamically equivalent to the billiard flow, we ask questions regarding the 
ergodic nature of the conjectured geodesic and billiard flows on the hypothesized 
'fractal flat surface' and the corresponding fractal billiard. We then state an open 
question asking whether or not it is possible to construct analogs of various classical 
trace formulas (e.g., the Gutzwiller trace formula) which can be used, in the classical 
billiard case, as a tool for connecting the length spectrum of particular billiards and 
the eigenvalue spectrum of the Laplacian defined on their associated fractal drums. 

In the long-term, we hope that the present (preliminary) study of the Koch 
snowflake billiard will help lay the foundations for a general theory of fractal bil- 
liards. 

1.1. Index of notation. At this point, we would like to provide the reader 
with an index of notation. Such an index will list the notation, a brief explanation 
of the notation and where this notation is first used or defined. 



I, [0,1] 


The unit interval. By convention, / is the base of the equi- 






lateral triangle A = KSq, and is identified with [0, 1] 




KS n 


The nth Koch snowflake prefractal curve, n > 


® 


KS 


The Koch snowflake fractal curve 


® 


n(Ks n ) 


The nth Koch snowflake prefractal rational billiard 


® 


Cl(KS) 


The Koch snowflake fractal billiard 


m 


!b 


The billiard map associated with the billiard fl(B) 


m 


fn 


The billiard map associated with the billiard fl(KS n ) 


m 


S(R) 


The flat surface corresponding to the rational billiard il(R) 


lh 


V 


The ternary Cantor set 


□j 


S(KS n ) 


The nth Koch snowflake prefractal flat surface 


on 




A side of Sl(KS n ) 


on 


~0 
•^n 


An initial basepoint of an orbit of ft(KS n ) 


M 


™k n 


A basepoint of an orbit of fl(KS n ) 


M 




An orbit of Sl(KS n ) with initial condition (a;°,0°) 


M 


9n,k 


A ghost of a side s n , k of &{KS n ) (1 < k < 3 • 4") 


M 


C n ,k 


A cell of Q(KS n ) corresponding to a side s n -i,k of 


m\ 




n(i<s n ^) 






A piecewise Fagnano orbit of ft(KS n ) 






A specific collection of points in the unit interval I 


m\ 


0.uim 2 ... 3 


A base-3 expansion of a number in the unit interval / 


m 


^n,k; ^ n ,k 


The ternary Cantor set of a side s n ^k of Q(KS n )\ the non- 


m 




ternary points of ^ n ,fc 




l,c,r 


An alternate alphabet used to give a ternary representation 


S21 




of an element of I 
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2. Background 

2.1. Mathematical Billiards and the Billiard Map. Under ideal condi- 
tions, we know that a point mass having a perfect elastic collision with a C 1 surface 
(or curve) will reflect at an angle which is equal to the incoming angle, both mea- 
sured relative to the normal at the point of collision. 

Consider a compact region fl(B) in the plane with connected boundary B. 
Then, Cl(B) is called a planar billiard when B is smooth enough to allow the Law 
of Reflection to hold, off of a set of measure zero. When B is a nontrivial connected 
polygon in K 2 , Q(B) is called a polygonal billiard, and the collection of vertices 
forms a finite set, which is a set of zero measure (when we take our measure to 
be the Hausdorff measure or simply, the arc-length measure on f2 (£?)). For the 
reader's easy reference, we provide a formal definition of rational billiard below. 

Definition 2.1 (Rational polygon and rational billiard). If B is a nontrivial con- 
nected polygon such that for each interior angle Oi of B there are relatively prime 
integers pi > and qi > such that 6i = — 7r, then we call B a rational polygon 
and fl(B) a rational billiard. 
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Figure 4. A billiard ball traverses the interior of a billiard and 
collides with the boundary. The velocity vector is pointed outward 
at the point of collision. The presence of a well-defined tangent 
at this point provides for the existence of a normal to the tangent 
and the recovery of the Law of Reflection (i.e., Snell's Law). The 
resulting direction of flow is found by either reflecting the vector 
through the tangent or by reflecting the incidence vector through 
the normal and reversing the direction of the vector. We use the 
former method throughout this paper. A rigorous discussion of the 
Law of Reflection in this context is given in |Sm] . 



The Law of Reflection essentially amounts to reflecting the incoming vector 
through the normal vector and then reversing the direction. Instead, we adhere 
to the convention in the field of mathematical billiards according to which the 
vector describing the position and velocity of the billiard ball (which amounts to 
the position and angle, since we are assuming unit speed) be reflected in the tangent 
to the point of incidence. Then we can rigorously reformulate the Law of Reflection 
as follows: the vector describing the motion is the reflection of the incoming vector 
through the tangent at the point of collision. Moreover, we can identify these 
two vectors and form an equivalence class of vectors in the unit tangent bundle 
corresponding to the billiard table Q(B). (See Figure 2] and [Sm for a detailed 
discussion of this equivalence relation on the unit tangent bundle tt(B) x S 1 .) 

For the remainder of this section, we suppose that fl(B) is a polygonal billiard 
with connected boundary!! Denote by S 1 the unit circle, which we may consider 
to represent all the possible directions (or angles) in which a billiard ball may 
initially move. In practice, one restricts one's attention to (B x <~ when 

discussing the phase space of the billiard fl(B). To clearly understand how one 
forms equivalence classes from elements of B x S 1 , we let (x, 9), (y, 7) G B x S 1 and 
say that (x, 9) ~ (y, 7) if and only if x = y and one of the following is true: 

(1) x — y is not a vertex of the boundary B and 9 = 7, 



In fact, in the remainder of the paper, except when B is the Koch snowfiake, B will be a 
nontrivial connected polygon. As such, the billiard flow given by the billiard map fg : (Q(B) X 
S 1 )/ ~ -> (Q(B) X S 1 )/ ~ will have finitely many singularities, whereas Q(KS) has, a priori, 
uncountably many singularities. 
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(2) x = y is not a vertex of the boundary B, but x — y is a point on a 
segment Si of the polygon B and 9 — 71(7), where r% denotes reflection in 
the segment Sj, 

(3) If x — y is a vertex of B, then we identify [x, 9) with (y,g(^f)) for every 
g in the group generated by reflections in the two adjacent sides having x 
(or y) as a common vertex. 

The billiard map (sometimes also called the Poincare map) is denoted by fg and 
is determined in the following way. Let [(x°,9 )} G (B x S 1 )/ ^ be the equivalence 
class of (x°,9°), with x° £ B and 9° e S 11 an inward pointing direction (unit 
vector) based at x°. Then /b([x o ,0 ]) identifies the first point of collision in the 
boundary and the outward pointing vector 9° that is parallel with the unit vector 
in the direction 9°. Since fg ■ (B x S 1 )/ ~ — > (B x S 1 )/ ~, we have that 
(a; 1 ,^ ) e [(a; 1 ,*? 1 )] = fg([x°, 6*°]). As expected, the outward pointing vector is 
identified with the inward pointing vector and one then takes as the direction of 
travel the inward pointing vector 9 1 based at x 1 ^ 

The billiard map fg describes the discrete billiard flow on the phase space 
(fi(B) x S 1 )/ ~. It is related as follows to the continuous billiard flow fg-.Rx 
(fi(B) x S 1 )/ ~ -> (Q(B) x S" 1 )/ ~. When we restrict our attention to B x S 1 , the 
flow f g on (B x S 1 )/ ~ is discrete. Specifically, one denotes the time at which the 
billiard ball collides with the boundary B by i 6 Z. The discrete billiard flow 
is then given by fg : Z x (B x S 1 )/ ~ — > (B x S* 1 )/ ~. Such a flow is also called 
the Poincare section of the continuous flow f B ; see |KaHal] for a more complete 
discussion of the billiard flow. In the sequel, instead of denoting the discrete flow 
by fg, we use the standard notation f B . In general, for every integer k > 0, we 
have [x k , 9 k ] = fs([x°, 0°]), where fg (defined above) can also be viewed as the fcth 
iterate of the billiard map fg . 

Remark 2.2. So as not to introduce cumbersome notation, when we begin dis- 
cussing the billiard map fj(S n corresponding to the nth prefractal billiard il(KS n ), 
we will simply write fxs n as /„. Moreover, when discussing the billiard flow on 
(tt(KS n ) x S 1 )/ ~, the fcth point in an orbit [{x k ,9 k )] € (Cl(KS n ) x S 1 )/ - will 
instead be denoted by [(a;*", 9 k ")], so as to be clear as to which space such a point 
belongs. 

In the event that a basepoint x 1 of fB{x°,9°) is a corner of 0(B) (that is, a 
vertex of the polygonal boundary B) and 9° was a direction for which the billiard 
flow would be periodic, then the resulting closed orbit is said to be singular. In 
addition, since 9° is a direction for which the resulting orbit is periodic, there exists 
a positive integer k such that the basepoint x~ k of fg k (x°, 9°) is a corner of f2(B). 
(Here, f B k denotes the fcth inverse iterate of fg.) The path then traced out by the 
billiard ball connecting x 1 and x~~ k is called a saddle connection. 

Remark 2.3. Within the subject of mathematical billiards, there appears to be a 
slight abuse of language. One may refer to the orbit of a billiard as the path traced 
out by the billiard ball or as the collection of incidence points in the boundary. In 
the latter case, such a set of points is referred to as the Poincare section of the 
billiard map fg- When we want to be clear as to which concept we are referring, 



The representative element is always the one for which the vector describing the direction 
after collision is inward pointing. 
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Figure 5. These two orbits are equivalent. The initial angle is 7r/3 
and the segments in the path remain parallel after each subsequent 
reflection. It should be noted, however, that an orbit with an initial 
basepoint at the midpoint of the base of the equilateral triangle and 
an initial angle of 7r/3 results in what is called the Fagnano orbit, 
which has period 3 and a path with a length that is half of the 
length of either orbit shown in this figure. This does not mean 
that we are making an exception and saying that this Fagnano 
orbit should be equivalent to the other two shown, but merely 
that one must be careful when making the claim that all orbits 
with the same direction are necessarily equivalent. 

we will specifically write 'the path corresponding to the orbit' or 'the footprint (or 
Poincare section) of the corresponding orbit', respectively. 

When considering periodic orbits (which are, by definition, non-singular closed 
orbits) , one point of interest is the corresponding path traced out by the billiard ball 
and the length of such a path. Such a path necessarily has finite length since the 
orbit has finite period. If two periodic orbits (?{x, 9) and &(y, cj>) have corresponding 
paths @*{x, 9) and &{y, 4>) that have equal lengths, equal periods and, x and y lie 
on the same segment of B such that &>(x, 9) and &(y, <fi) remain parallel from one 
basepoint to the next, then we say that these two orbits are equivalent. From this, 
we see that equivalence of orbits does not depend at all on the choice of the initial 
basepoint of the given periodic orbitsQ see Figure [5J 

2.2. Flat Structures and Flat Surfaces. In this section, we deal only with 
flat surfaces constructed from rational billiards; see Definition 12. II in §2.1l above. By 
showing that there is an intimate connection between what is called a flat surface 
and a billiard Sl(B), we will eventually be able to demonstrate that reflection in 
certain corners can be defined and, for certain rational billiards, the nature of an 
orbit of a special class of rational billiards does not depend at all on the initial 
point, but only on the initial direction^ To such end, we discuss the notion of a 
flat structure. 



^Except for a few special cases, one of which is discussed in the caption of Figure [5] 
In Jj6j we discuss what are called uniquely ergodic flows and the Veech dichotomy. Un- 
derstanding such concepts is not necessary for §§3-5; so we defer explanation of those concepts 
until §^6]and refer the reader to [Vel— 3] for an immediate and detailed description of the Veech 
dichotomy and to MasTa, HuSc, Zo for a pedagogical discussion of uniquely ergodic, periodic 
and aperiodic billiard flows. 
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A flat structure on a connected surface M is an atlas (U a ,4>oi) with a finite 
number of singularities such that, away from these singularities, each 'coordinate 
changing' map 

(2.1) ^o^ 1 ■. <j)f3(u a nup) ^ <jy a (u a nu fj ) 

is a translation in R 2 . Specifically, a flat structure satisfies the following definition. 

Definition 2.4 (Flat structure). Let M be a compact, connected, orientable sur- 
face. A flat structure on M is an atlas u, consisting of charts of the form (U a , <f> a ), 
where U a is a domain (i.e., a connected open set) in M and <p a is a homeomorphism 
from U a to a domain in R 2 , such that the following conditions hold: 

(1) the domains U a cover the whole surface M except for finitely many points 
Z\, Z2, Zk, called singular points] 

(2) all coordinate changing functions are shifts (i.e., translations) in R 2 ; 

(3) the atlas u> is maximal with respect to properties (1) and (2); 

(4) for each singular point Zi, there is a positive integer rm, a punctured 
neighborhood Ui not containing other singular points, and a map ipi from 
this neighborhood to a punctured neighborhood Vi of a point in R 2 that 
is a shift in the local coordinates from co, and is such that each point in 
Vi has exactly m, preimages under 

Definition 2.5 (Flat surface). We say that a connected, compact surface equipped 
with a flat structure is a flat surface. 

Remark 2.6. Calling a connected, compact (2-dimensional) manifold with flat 
structure a flat surface is somewhat of an abuse of language, but it enables us to 
be briefer and refer to related notions with greater ease. Note that in the literature 
on billiards and dynamical systems, the terminology and definitions pertaining to 
this topic are not completely uniform; see, for example, [GaStVo,Gul,GuJul 
2,HuSc,KaHa2,Mas,MasTa, Vel-3,Vo,Zo]. We have adopted the above defi- 
nition for clarity and the reader's convenience. 

The singular points of a flat surface S are called singularities of the flow. There 
are two types of singularities in a flat surface: removable and nonremovable. They 
are called such because it may or may not be possible to define the flow at these 
points. We are primarily interested in how the geodesic flow behaves at these points. 
We now turn to a discussion of why these singularities can be termed "removable" 
or "nonremovable" and how to discern between the two types. Consider Zi in the 
set of singularities of a flat structure. This point has what is called a conic angle. 
A conic angle is the number of radians required to form a closed circle about 

We state the following definitions, for completeness. 

Definition 2.7 (Removable conic singularity). A singularity of a flat structure 
(U a , <p a ) is a removable singularity of the flow if the conic angle about such a point 
is 2-7T. 

Definition 2.8 (Nonremovable conic singularity). A singularity of a flat structure 
(U a ,<ft a ) is a nonremovable singularity of the flow if the conic angle about such a 
point is 27rc, for some integer c > 1. 
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With regards to the geodesic flow on a flat surface S, removable singularities 
of the flow pose no problem and the flow lines through such points can continue 
unimpeded. However, when a singularity is nonremovable, the flow cannot continue 
unimpeded. Moreover, when zi is a removable singularity of a flat structure uj, there 
exists a flat structure Co on S such that Zi is not a singular point of Co. Technically 
speaking, when Zi is a removable singularity of the flow on a surface with a flat 
structure, there exist neighborhoods Up and Uj of Zj and maps <fip an< i fty such 
that the associated transition map <f>p o tfi^ 1 is a shift in the local coordinates at 
Zi. Consequently, the geodesic flow may be continuously extended at removable 
singularities of the flat surface. 

We now discuss how to construct a flat surface from a rational billiard. Consider 
a rational polygon billiard fi(P) with m sides and interior angles ir(ui/vi) at each 
vertex z%, for 1 < i < m. Here, m and are relatively prime positive integers. 
Then, a (laborious) calculation shows that for some j < m, Vj = lcm(wi)™ =1 f jy. 
Consequently, the linear portions of the planar symmetries generated by reflection 
in the sides of the polygonal billiard fi(P) generate the dihedral group Dn, where 
N := Icmlvi}™^. Here, by definition, the dihedral group Dn denotes the group of 
symmetries of the regular iV-gon. So as to be clear, we mention that the notation 
Dn does not refer to the (wrong) fact that such a finite group has N elementsQ 
We next consider f2(P) x £)jv (equipped with the product topology). We want to 
glue 'sides' of f2(P) x Dn together and construct a natural atlas on the resulting 
surface M so that M becomes a flat surface. 

To such end, let p\ — p2 be a point on a side s a of H(P), r a be the linear portion 
of the reflection determined by reflecting fi(P) in the side a, and {p\, r%), (p 2 , r 2 ) £ 
P x Dn- Then, by definition, (pi, ri) ~ (p2, r 2 ) if and only if 

(1) (pi,ri) = (p 2 ,r 2 ), or 

(2) r a = rf V 2 , or 

(3) pi and p 2 are the same vertex of f2(P) having adjacent sides s a and Sb, 
with r^ x r 2 belonging to the subgroup generated by r a and rb- 

It is easy to check that ~ is an equivalence relation on f2(P) x Dn. More work 
is required to show that M := (O(P) x Dn)/ ~ is a compact, connected, orientable 
surface. As a result of the identification, the points of M that correspond to the 
vertices of f2(P) constitute (removable or nonremovable) conic singularities of this 
surface. Heuristically, f2(P) x Dn can be represented as {rjP}?^, in which case it 
is easy to see what points are made equivalent under the action of ~. 

Denote by f2(P)° the interior of the billiard fi(P). To construct a flat structure 
on M, let L/i = ri(P) x {n}. Then {[/i,0j}i=i can be naturally extended to 
constitute a flat structure on M, in the sense of Definition 12.41 Hence, M is a flat 
surface, in the sense of Definition ^. 51 The map <j) i o(j)'J 1 : cf>j (Ui CiUj) — !• cf>i{Ui fl C/j) 
is a translation in the local coordinates of a point z € UidUj ; i.e., ipiO^J 1 ^) = z+d, 
where the constant d is independent of the choice of i and j'0 



^Actually, Dn has 2JV elements, and the standard group theory notation for the dihedral 
group is then -D2JV1 since the cardinality of the group is often given more importance, from the 
perspective of group theory. 

^A priori, the choice of d describing the translation in the local coordinates does depend on 
i and j. However, given the fact that one constructs the flat surface by identifying parallel and 
opposite sides of the polygon B, for a fixed direction 8, one can describe a parallel line field in the 
direction 8 := arctan^-, with d =: a + \/— 1/3. 
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Example 2.9 (The equilateral triangle billiard Cl(KSo)). We consider the equilat- 
eral triangle A := KSo, an important example in the literature (see, e.g., [BaxUmj ) 
and an even more important example when considering the Koch snowflake prefrac- 
tal billiard. Since the interior angles Ttpi/qi are all the same and equal to ir/3, we 
have 



(2.2) N = lcm(gi)? =1 = lcm(3, 3, 3) = 3. 

Consequently, the associated surface is given by S(KSo) := (AxD^)/ ~. Moreover, 
there is a flat structure on this surface and since all the singularities are removable 
conic singularities, such a structure can be extended to include the singular points 
and the resulting surface is topologically equivalent to a torus. (Really, the associ- 
ated flat surface S(K Sq) is the hexagonal torus, which is topologically equivalent to 
the standard torus; see Figured in §5. ) 

We close this discussion by recalling an important fact about the geodesic flow 
on a flat surface S(B) constructed from a rational billiard il(B) and the billiard 
flow on that rational billiard. These two flows can be shown to be (dynamically) 
equivalent under the action of the group Dn associated with the construction of the 
corresponding flat surface. Heuristically, one may view the corresponding equiva- 
lence as follows: a billiard flow line may be straightened to a geodesic flow line and 
additionally, a geodesic flow line may be collapsed into a billiard flow line. In order 
to be more technically correct, we further explain the details of the equivalence of 
these two flows. If we consider the geodesic flow on the surface, the quotient of 
the phase space by the group of symmetries associated with the construction of 
S(B) has the effect of collapsing the space to a space with a quotient flow that is 
isomorphic to the billiard flow. On the other hand, any given billiard orbit may be 
straightened by making successive reflections, via the action of Dm, in the identi- 
fied sides of the flat surface, therefore producing a straight-line flow line on the flat 
surface S{B). 

2.3. Inverse limit sequence and inverse limit. Inverse limits of various 
topological (algebraic, or geometric) objects will play an important role in this 
paper. Hence, since such a notion may not be familiar to all readers, it may be 
helpful to recall some basic facts pertaining to this subject. Further information can 
be found, for example, in |HoYoj and, in a more general context, in |Bo[ IMcLj . 

We discuss the inverse limit in the context of the category Set. The objects 
are sets and the morphisms are set maps@ Consider a partially ordered collection 
of sets {Xi}^ 1: with each Xi equipped (for each integer n > 1) with a map F n : 
n°^i Xi — > X n defined by F n ((xi)^.i) = x n . If for every n and every m < n, where 
m,n € N* := {1,2,...}, we define the map F nm : X n -> X rn by F nrn (x n ) = x m , 
then the set 



(2.3) UmXj := I (x,)^ £ X t \F nrn (x n ) = x m , for all m < n 



"Discussing the inverse limit in the context of Set is purely a formality that allows us to 
speak in more concrete terms and utilize important existences and uniqueness properties of the 
inverse limit in the category Set. See |Bo[ IMcLI for a general discussion. 
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exists and is unique; it is called the inverse limit of the inverse limit sequence 
{Xi,Fi\^ l . Furthermore, the maps F nm are called the transition maps of the 
inverse limit system. 

Naturally, if we work instead within the category of Topological Spaces, then 
all the maps involved should be morphisms of that same category (and hence, here, 
continuous maps). 

We next give an example of the inverse limit of a particular collection of sets 
(or rather, of topological spaces). 

Example 2.10 (The ternary Cantor set %?). Recall that a ternary number is a 
number expressed in terms of a base-3 number system, say, in terms of the char- 
acters (or symbols) {0,1, 2}0 For example, 1102 in base-3 is actually the number 
38 in the base-10 number system. When one constructs the Cantor set ^ , one ma 
do so by removing middle thirds ( open intervals ) of successive approximations)^ 
Let n > 1. In removing the middle third from an interval of length 1/3™ _1 ; we are 
essentially producing a left third and a right third interval of length 1/3™. As such, 
we can label the left third as and the right third as 2. One quickly sees that the 
Cantor set is a collection of infinite words written entirely in terms of 0's and 
2's. There are no ternary numbers in ^ whose address contains 1, because that 
would mean that we did not remove a middle third from some interval in the con- 
struction processV^ More importantly, our labeling system described above enables 
us to determine particular ternary numbers. 

Let be the nth prefractal approximation of the Cantor set. In the context 
of the inverse limit construction and addressing system above, is the collection 
of 2™ points, each having an address given by a finite ternary expansion of length 
n and each never containing the character 1. On the other hand, in the context 
of the geometric construction detailed in the previous paragraph, ^„ consists of 2™ 
compact intervals of length 1/3™ (i.e., of 'scale n'). Furthermore, each address 
in %? n may be thought of as an address of a particular segment that remains after 
removing 2™ _1 segments from the unit interval I; see Fiaure lUt^ We take the more 
geometric interpretation as the definition of^ n ; in that case, the sequence {^n}^! 
of prefractal approximations converges to the Cantor set ^ : %f„ — > c € as n — > oo. 
in the sense of the Hausdorff metric. (Also, more simply, c (o n is monotonically 
decreasing and c € = H^Li ^n-) We will next discuss another way in which ^ can 
be viewed as the 'limit' of { < Tg n }'^L 1 . 



^Any three symbols suffice. One can, and we do so in J5] represent elements of the Cantor 
set in terms of the characters {i,c, r}, where such characters stand for left, center and right, 
respectively. In general, any element in the unit interval / can be represented by a finite or an 
infinite sequence expressed in terms of an alphabet consisting of the characters l,c,r. 

^This is not the only way to construct the ternary Cantor set, but most, if not all, methods 
amount essentially to the same process. 

^For the endpoints of the deleted intervals (also called 'ternary points' and necessarily of the 
form p/3 q , for p, q nonnegativc integers with p < 3 q and p not divisible by 3, when we restrict our 
attention to the unit interval I), the resulting address is not unique and may contain l's, although 
only finitely many l's. For example, one may represent 1/3 by the finite ternary expansion 
0.1 or by the infinite ternary expansion 0.02, where the overbar indicates that 2 is repeated ad 
infinitum. We adhere to the convention that such ternary numbers are always represented by 
infinite expansions given in terms of only 0's and 2's. 

12 Such a construction is called construction by tremas, where one removes segments ad 
infinitum, thereby producing the fractal set. At each stage, middle thirds are removed from the 
remaining segments, thereby producing the Cantor set in the limit. 
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Figure 6. The geometric construction of the Cantor set. One 
removes middle thirds from remaining segments until all that re- 
mains are the points of the Cantor set. Benoit Mandelbrot is cred- 
ited with the coining of the illustrative name Cantor Dust, a name 
that is quite fitting as it appears that all that remains in the limit 
is 'dust.' Shown here in this figure are the first six approximations 
of the Cantor set. 

If for positive integers m < n, we now define the transition map T nm : c <f n — > 
^ m as the truncation map that truncates addresses of segments in the prefractal 
approximations c £ n to addresses of segments in the prefractal approximation %f m 
by simply removing the last n — m characters from the address, then we form an 
inverse limit of the prefractal approximations ^ n that is exactly the Cantor set < H> : 



More background information about inverse limits is provided, for example, 
in [HoYol §2-14 and §2-15], where the following well-known theorems can be 
found; sec [H o Yo) . Theorems 2-95 and 2-97, along with Corollaries 2-98 and 2- 
99. (All the topological spaces considered in Theorems 12.111 12.121 and 12.131 below 
are implicitly assumed to be metrizable. Moreover, by the "Cantor set" , we mean 
the classic ternary Cantor set discussed in Example 12.101 just above.) 

Theorem 2.11. The inverse limit of finite sets is a compact and totally discon- 
nected set. Conversely, any such topological space is homeomorhpic to an inverse 
limit of finite sets. 

Theorem 2.12. Any compact and totally disconnected space is homeomorphic to 
a (closed) subset of the Cantor set. 

Theorem 2.13. Any two totally disconnected and verfec¥^ compact spaces are 
homeomorphic to one another (and hence also to the Cantor set). 

Remark 2.14. In the literature on dynamical systems, it is common to use the 
term topological Cantor set to refer to a totally disconnected and perfect compact 
space (i.e., to a metrizable space that is homeomorphic to the Cantor set). 

Remark 2.15. In §5, we will show that what we will be referring to as the foot- 
prints of the primary piecewise Fagnano orbit and, in general, piecewise Fagnano 

1 ^Recall that a subset of a topological space is called perfect if it is closed and contains no 
isolated points. Hence, a compact space is perfect if it has no isolated points or equivalently, if 
each of its points is a limit point. 




c m for all m < n 
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orbits of the Koch snowflake billiard il(KS) are, in fact, topological Cantor sets 
(see Thcorcm l5.5[) . Moreover, the footprint of the primary piecewise Fagnano orbit 
will be the analog of the classic ternary Cantor set, and as a subset of the Koch 
snowflake KS, the union of the footprints of every piecewise Fagnano orbit will con- 
stitute a subset of what we will call the elusive limit points of the Koch snowflake 
KS. 

3. The Flat Surface S(KS n ) as a Branched Cover of S(KS ) 

We denote the flat surface M as constructed from a particular rational billiard 
f2(P) by S(P). In particular, S(KS n ) is the flat surface associated with the pre- 
fractal billiard il(KS n ). The flat surfaces S(KS n ), n — 1,2, 3, are given in Figure 
[7] For each billiard ft(KS n ), the group of symmetries Djv, where N — lcmjwj}^ 
(that is, the second component in the product Q,(KS n ) x Dm) is the dyhedral 
group -D3, and thus is independent of n. From this, we deduce that for any 
n > 0, there are six copies of the prefractal billiard table ft(KS n ) (with sides 
appropriately identified) used in the construction of the associated flat surface 
S(KS n ) := (Cl(KS n ) x D 3 )/ ~; see Figure [7] We refer the reader back to 
for the discussion of flat surfaces and the associated conical singularities. 

Definition 3.1 (Covering map). Let E and B be topological spaces. A covering 
map p : E — > B is a continuous surjective map such that each point b € B admits 
a neighborhood U of b for which the preimage j> _1 (f7) is a disjoint collection of 
open sets in E, each of which is mapped homeomorphically onto U via p. One then 
says that U is evenly covered by p, and that E is a covering space for B; see, for 
example, Ma, Chap. 5]. 

Definition 3.2 (Branched (or ramified) cover). Let E and B be topological 
spaces. A continuous map p : E — > B is a branched cover of B if for all but a 
finite number of points of B, p is a covering map of E onto B. The set of points of 
B that are not evenly covered by p is called the branch locus (or set of ramification 
points). 

Example 3.3 (The map p : C — !• C). The map p : C — > C, given by p(z) = z 2 , is a 
branched covering ofC, with branch locus {0}. Hence, it is certainly not a covering 
map. On the other hand, p : C — {0} — > C — {0}, given by the same expression 
p{z) — z 2 , is a covering map, since it is locally trivial: indeed, each nonzero complex 
number z in the target space has an open neighborhood U such that p, restricted to 
p _1 (J7), is equivalent to the projection onto U x {+, — }. 

For the remainder of the paper, when we say that a regular polygon is of scale 
n, we mean that the side length of the regular polygon is 1/3™. For example, an 
equilateral triangle of scale n is one for which the side length is 1/3™. 

Taking as inspiration the results and methods of Gutkin and Judge in [GuJul] 
and [GuJu2j . we now show that for each n > 1, the flat surface S(KS n ) is a 
branched cover of the hexagonal torus S(KSo); see Corollarv l3.7l To such end, we 
establish several results culminating in the fact that S(KS n ) is tiled by equilateral 
triangles of scale n. 

Lemma 3.4. Let n e N. Then, for any positive integer k > n, S(KS n ) can be 
tiled by equilateral triangles of scale k. 
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Figure 7. The flat surfaces S(KS l ), i = 1,2,3. Note that the 
proper identification is not shown in the figures above. Given the 
arrangement of the six copies of KS n , one then identifies opposite 
and parallel sides to make the proper identification that results in 
a geodesic flow that is dynamically equivalent with the billiard flow 
on the associated billiards fl(KSi),fl(KS 2 )M(KS 3 ). 
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Figure 8. We see that A„ is tiled by nine copies of A n +i, with a 
hexagonal tile H n+ \ in the center. 




Figure 9. The hexagonal torus S(KSq). (As usual, similarly 
marked sides are identified.) It should be noted that S(KSq) is 
topologically (but not metrically) equivalent to the flat square 
torus. 

Proof. This follows from the construction of the Koch snowflake. We note 
that each triangle of scale n, denoted A„, can be tiled by Q k ~ n triangles of scale k > 
n; see Figure[5]for the case when k = n + 1. Note that S(KS n ) — (fl(KS n ) x D 3 )/ ~ 
and that £l(KS n ) is constructed from £l(KS n ~i) by gluing a copy of the billiard 
table fl(KS n ) of scale n to every side Sk at the middle third of Sk and then removing 
the segment common to Sk and A n+ i. So, S(KS n ) can be tiled by equilateral 
triangles of scale k. □ 

In the sequel, given a bounded set A C R 2 , we will write that 11 A can be tiled by 
H n " in order to indicate that A can be tiled by finitely many copies of hexagonal 
tiles H n of scale n. 

Lemma 3.5. Let n G N. Then the hexagonal torus S(KSq) can be tiled by H n . 

Proof. Consider the hexagonal torus S(KSq), as shown in Figure|Hl In Figure 
[TOl we see that S(KSo) can be tiled by nine hexagons.. As mentioned in the proof 
of Lemma 13.41 A„ can be tiled by nine copies of A„ + i overlapping only at the 
edges. Moreover, at the center of each A„ is a hexagon of scale n + 1, H n+ i; see 
Figure [U 

A hexagon H n contains six equilateral triangles A„. As we have seen before, 
each triangle A„ has at its center a hexagon H n+ i; see FigureEJ These six triangles 
A„ are then all tiled so that for each A„, there is a A n+ i such that a total of six 
copies of A„_|_i comprise an additional hexagonal tile H n+ i placed at the center of 
A„. Consequently, H n is tiled as shown in Figure [TT] Then, as shown in Figure IT2| 
three copies of H n arrange so that at each common vertex, there are three rhombic 
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Figure 10. We see that there are seven hexagonal tiles of scale 
n = 1 tiling S(KSq). By means of the identification, we recover 
the other two hexagonal tiles, as claimed. 




Figure 11. Six triangles A„ tile H n . The hexagonal tile H n is 
tiled by seven tiles H n+ i with six rhombic tiles. 

tiles comprising another additional copy of H n +i- Therefore, for any n, S(KSq) is 
tiled by hexagons H n . □ 

Proposition 3.6. For every neN, S{KS n ) can be tiled by H n+ i in such a way 
that each conical singularity is at the center of a hexagonal tile. 

PROOF. We proceed by induction on n. First, assume that n = 1. Then 
S(KSi) can be tiled by H 2 so that every conic singularity of S(KS\) is at the 
center of some hexagonal tile. Now, given n > 2, suppose that S(KS n -x) is tiled 
by hexagons of scale n so that every conic singularity is at the center of some 
hexagonal tile of scale n. 

We first note that the nature of a conic singularity in the surface dictates that 
such a singularity is actually common to four hexagonal tiles. We now embed 
in S(KS n ) the tiling of S(KS n -i) by H n . We see by way of the equivalence 
relation on Q(KS n ) x D n that for every triangular region A„ in S{KS n ) that is 
not in S{K SVi-i), there are five other regions in S{KS n ) that comprise a complete 
hexagonal tile H n , such that every side of the hexagonal tile is common to some 
other hexagonal tile in the tiling embedded from S(K S n -i). Furthermore, by 
Lemma 13.51 each hexagon can be tiled as shown in Figure 1111 

Note that every hexagon H n is bordered by six other hexagons H n . Therefore, 
by the proof of Lemma l3.4( every vertex of H n is at the center of a hexagon H n+ \. 
Moreover, vertices of H n with edges common to a segment of the polygon KS n 
coincide with ternary points ^ft, where p = 1 (mod 3) or p — 2 (mod 3). 
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Figure 12. We see that three hexagonal tiles H n tiled as in Figure 
1111 can be arranged so that a rhombic tile from each all combine 
to form another hexagonal tile -ffn+i- 



Consequently, vertices of any tile H n+ i with a side common to £l(KS n ) that 
are not at a conic singularity of S(KS n ) are a distance 3 X+i from a conic singularity 
of S(KS n ). This is exactly the distance to the center of a hexagon tile H n+ i along 
side of A„ + i tiling H n+ \. Hence, we conclude that every conic singularity is at the 
center of some hexagonal tile H n +i, as desired. □ 

Corollary 3.7. For every n€N, the prefractal Koch snowflake flat surface S(KS n ) 
is a branched cover of the prefractal Koch snowflake flat surface S(KSq). 

Proof. The center point Xq of the flat hexagonal torus S(KSq) is a branched 
locus of the cover S(KS n ) when S(KS n ) is tiled by H n+ i as described in Proposi- 
tion 13.61 This follows from the fact that every conic singularity is at the center of 
four hexagonal tiles. Specifically, this means that this center point zn is not evenly 
covered by the quotient map p n : S(KS n ) — > S(K S n ) / (3 n+1 H n+ i)\^j Any other 
point in S(KSo) is evenly covered since every element in the fiber p~ 1 (z), z ^ zq, 
has a conic angle of 2ir. □ 

We close this section by highlighting a possible pattern in the construction of the 
flat surfaces S(KS n ). In Figure UM we see that, under the proper re-identification, 
three equilateral triangle tiles appended to each copy of the equilateral triangle 
in the flat surface S(KSq) (see Figure [T3"|) results in the flat surface S(KSi). In 
Figure 1151 the proper re-identification yields three tori that are interconnected in 
such a way that results in a flat surface with genus g = 10. Heuristically, one tears 
three holes in the hexagonal torus, then glues three additional tori to the existing 
flat surface in such a way that the proper surface results. While this is admittedly 
very difficult to visualize, one can get the impression that each surface S(KS n ) 



^The notation 3 n+1 -ff n +i indicates that we are scaling the hexagonal torus of scale n + 1 
by 3 n+1 , thereby producing Hq, the hexagonal torus. 
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Figure 13. Six equilateral triangles glued together correctly con- 
stitute a hexagonal torus. We show the surface in an exploded view 
so as to emphasize the fact that there are six copies of the equilat- 
eral triangle embedded in the surface. Prior to identifying properly, 
this figure constitutes the heuristic description of Q(KSq) x D3. 




Figure 14. Here we show the surface S(KSi) in an exploded 
view. Opposite sides that are parallel and opposite are identified. 
We also see that the surface S(KSi) can be viewed as gluing three 
tori to the existing torus in such a way that the resulting flat surface 
has genus g — 10. 



results from S(KS n —i) by tearing and gluing to the newly opened holes 3 • 4™ many 
equilateral triangles of scale n to S(KS n -i)- 
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Figure 15. We show here the three tori that are appended to the 
existing torus, rearranged to clearly see the three tori discussed in 
the caption of Figure HU 

4. Periodic orbits of f2(KS n ) in the direction of ir/3 

In the previous section, we saw that S(KS n ) is a branched cover of S(KSq). 
This implies that for every n £ N, a closed geodesic on S(KS n ) will project down to 
a closed geodesic on the hexagonal torus S(KSq) under the action of the covering 
map p n defined in the proof of Corollary 13. 71 Also, a closed geodesic 7 on S(KSq) 
lifts to a segment on S(KS n ) (that is not necessarily closed). However, according 
to the discussion in Sj3j there exists a positive integer k such that the lift of 7 fc is 
a closed geodesic on 6>(ir>SVi)I3 Since the geodesic flow on S(KS n ) is dynamically 
equivalent to the billiard flow on (il(KS n ) x S 1 )/ ~, it follows that a direction 
giving rise to a closed orbit in Vt(KSo) is a direction giving rise to a closed orbit 
in Q(KS n ) for every n > 0, and vice-versa. We may put this more succinctly as 
periodic directions in S(KSq) are exactly the periodic directions in S(KS n ), and 
vice-versa. We summarize the above discussion in the following theorem. 

Theorem 4.1. The geodesic flow on S(KSq) is closed if and only if for every 
n > 0, the geodesic flow on S(KS n ) is closed. Moreover, the set of directions for 
which the discrete billiard flow f n on VL{KS n ) is closed ^ {i.e., regardless of the 
initial basepoint, a direction for which a geodesic will be closed) is exactly the set 
of directions for which the billiard flow fo is closed on fl(KSo). 

Remark 4.2. If {ei, e-i\ is a basis for R 2 , then a vector z £ R 2 is called rational 
with respect to {ei,e2} if z = ne\ + me2, n/m € Q (that is, n,m G Z, m 7^ 
0). The plane can be tiled by VL(KSq). As proved in |Gu2| . the collection of 
directions that gives rise to closed orbits of the equilateral triangle billiard ft(KSo) 
is exactly the set of directions that are rational with respect to the basis {ei, e-i\ — 
{(l,0),(l/2, v / 3/2)}. By Theorem SU for every n > 0, the same collection of 
rational directions (with respect to {ei,e2}) describes the directions for which the 
billiard flow on il(KS n ) is closed. 

In the sequel, when we say that an angle is measured with respect to a fixed 
coordinate system, we mean that every direction is measured with respect to the 
same coordinate system up to translation, but not rotation. For example, a direc- 
tion of 7r/3 measured with respect to a side of a prefractal billiard table may not 
be 7r/3 when measured with respect to some fixed coordinate system. Rather, it 
may be 57r/3; see Figure [TBI In order to maintain a consistent measurement, we 



1 ^Recall that the notation 7* is meant to represent k — 1 many concatenations of 7 with 
itself: 7*7*...*7 = 7 fe 

^tt should be noted that we are making a slight abuse of notation and language. The billiard 
flow is a flow on the phase space (Q(KS n ) X S 1 )/ ~ and f n is the billiard map. Iterates of the 
billiard map then yield elements of the Poincare section. When we say the billiard flow is closed, 
we mean that the Poincare section is finite and vice-versa. 
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Figure 16. While each of these angles can be measured as 7r/3 
or 27r/3, relative to the sides from which each one emanates, we 
instead measure relative to a fixed coordinate system given in Fig- 
ure I17[ resulting in the measurement of each angle being equal to 
5tt/3. 

fix a coordinate system at the base of the equilateral triangle billiard fl(KSo); see 
Figure IT71 

Remark 4.3. Let n > and k < 3 • 4™. We denote by s n ,k a side of the billiard 
Q(KS n ). A basepoint of an iterate of the billiard map /„ := fKS n of Q(KS n ) is 
then denoted by x^ n . So as to be perfectly clear, we mention that the superscript 
k n is not related to the number of sides of fl(KS n ). Rather, it is a notation that 
helps us differentiate between 'fc' iterates of the billiard map /„ and / m , m ^ n. 
Hence, /*" is the k n th iterate of the billiard map /„. An initial basepoint of an 
orbit is then denoted by x° and an initial direction is 

Consider an inward pointing vector with basepoint x^ and angle 0" measured 
relative to the side of £l(KS n ) on which x„ resides (assuming x„ is not a vertex). 
Since we always want to measure angles relative to a fixed coordinate system, we 
denote the measure of the angle 0" relative to the fixed coordinate system by #(</>) ° . 
Because of the unique geometry of the prefractal approximation, for every inward 
pointing vector forming an angle <^>„ with a side s n ^ on which the basepoint lies, 
there exists k G Z such that @(<fi)n — kir/3± </>Fl 

We have the following definitions and theorems (see Figure[T5]for an illustration 
of Definitions Q3] and I4.5[) : 

Definition 4.4 (Compatible sequence of orbits). Let i^o(^O'^o) be an orbit of 
the equilateral triangle billiard ft(KSo) = f2(A). Assume that 0°) is the initial 
condition of an orbit £? n (x°,0°) of Sl(KS n ) such that 0° = 0§ (where each angle 
is measured relative to the fixed coordinate system) and a;° is collinear with Xq in 
the direction 9^ and there are no points of KS n between x^ and Xq. Then we say 
that ^„(a;°,0°) is compatible with ^o(^oi^o) anc ^ the resulting collection of orbits 
satisfying such a condition is called a compatible sequence of orbits. We write the 
compatible sequence as {^i(x°,9f)}°^ . 



''Hence, the angle 6(tt/3) is merely an integer multiple of 7r/3. 
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Figure 17. The fixed coordinate system relative to which all an- 
gles are measured. 




Figure 18. A compatible sequence of orbits. The first three or- 
bits in a compatible sequence are illustrated, and the first three 
initial basepoints in the compatible sequence of initial basepoints 
are indicated by the arrows. Here, the initial direction is 7r/3 and 
x° = 1/2. 



Definition 4.5 (Compatible sequence of initial conditions). Let 8® be a fixed 
direction measured relative to the fixed coordinate system. Assume that Xq is a 
point on the equilateral triangle KSq such that 8® based at x[j is inward pointing. 
Further, assume that x°, n > 1, is a point on the boundary of £l(KS n ) such that 
x° and Xq are collinear in the direction 8q with no points of KS n common to the 
segment (interior to 0,(KS n )) joining x[J and x°, when x$ is viewed as a point of 
il(K S n ) (if n — 0, then a;° is trivially collinear with Xq). Then we say that x° and 
Xq are compatible and that the sequence of compatible points {x°}°^ is a compatible 
sequence of initial basepoints; furthermore, if for each i > 0, 8® = 9q is an angle 
in S 1 (measured relative to the fixed coordinate system) that corresponds to an 
inward pointing vector at a basepoint x®, then {{x®,8®)}°^ Q is called a compatible 
sequence of initial conditions. 
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Theorem 4.6. A sequence of initial conditions {{x® , 0°)}£2. o is a compatible se- 
quence of initial conditions if and only if the corresponding sequence of orbits 
o is a compatible sequence of orbits. 

PROOF. This follows easily from Definitions 14.41 and [4.51 □ 

Theorem 4.7. If &q(xq,9q) is a closed orbit in Q,(KSq), then every orbit in the 
corresponding compatible sequence of orbits {&i(x®, 0°)}£2. o is a closed orbit. 

PROOF. Consider a closed orbit o (xq,8q) of ft(KSo). Then the period of 
^bC^oi^o) i s nn ite and the unfolding of the orbit in S(KSq) is a closed geodesic 
(or, in the case of the Fagnano orbit, twice the unfolding of the Fagnano orbit is a 
closed geodesic in S(KSq)). Then, by the discussion in §3, the lift of the geodesic 
to S(KS n ) is some segment 7 such that for some positive integer k, j k is a closed 
geodesic in S(KS n ). Then, the corresponding billiard orbit ^(x^, 0°) on Q(KS n ) 
is also closed (by the dynamical equivalence between the billiard flow on fl(KS n ) 
and the geodesic flow on the flat surface S(KS n )). 

The orbit & n (x° , 0°) is then compatible with &q(x% 0g). Since n was arbitrary, 
it follows that {^(x , 9^)}°^L is a compatible sequence of closed orbits. □ 

Definition 4.8 (Compatible sequences of closed and periodic orbits). If every 
orbit in a compatible sequence of orbits is closed, then we say that it is a compatible 
sequence of closed orbits. If, in addition, no orbit in the compatible sequence is 
singular, then we call the sequence a compatible sequence of periodic orbits. 

Remark 4.9. Let {^s(j/ , 0(0)°)}£2. o be a compatible sequence of orbits. If we 
want to discuss a particular compatible sequence of orbits, we use the fact that 
0(0)9 = 0(0)0 for all % > 1 and 1) write 0(0)°, as 0(0) and 2) ^(2/°, 0(0)°) as 
& n {Vni 0(0))- in the event that we are discussing a compatible sequence of orbits 
with 0(0)0 = 0(vr/3) for every n > 0, we write n (y° n , 0(0)°) as ff n [y° n , 0(tt/3)). 

Remark 4.10. The following definitions, Definitions l4.IIH4.i3l are also given 
in [LaNiell §3 & §4], up to a few minor differences which have allowed for the 
convenient definition of piecewise Fagnano orbit, which is given in Definition 14.141 

Definition 4.11 (Ghosts of KS n ). Let n > and {s n .fc}fc=i be the collection of 
segments comprising the polygonal boundary KS n of the billiard Q(KS n ). Then, 
for 1 < k < 3 • 4™, the open middle third of the side s Ui k is denoted by g n ^ and is 
called the ghost of the side s rh k- Moreover, the collection G n = {g n ,k}k=i is called 
the ghost set of KS n . The segments g n ^ are removed in order to generate KS n -\-i; 
see Figures H9Ta)-(c). 

Definition 4.12 (A cell C n> k of ft(KS n )). Consider (the 'set-theoretic differ- 
ence') £l(KS n ) \ Q(KS n -x). The resulting triangular regions are then called cells 
of il(KSn). We denote a cell of Q(KS n ) by C ni k, where k denotes the side of 
Sl(KSn-i) to which the cell was glued; see Figure 

Definition 4.13 (Ghost of a cell C n>k ). Let n > 1 and I < k < 3-4™- 1 . EC n , k is 
a cell of Q(KS n ), then the ghost g n -i,k corresponding to the side s n -i,k that was 
removed in the construction of Q(KS n ) is called the ghost of the cell C n ,k- 



18 Hence, there are 3 • 4™- 1 cells C„ ife of Q(KS n ) and so 1 < k < 3 ■ 4"" 1 , in both Definitions 
l4~12l and l4~T3l 
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A 



A 



(a) The ghost set of KS Q 
denoted by Go- 



= A, 



(b) The elements of the ghost set 
Go are removed. 




(c) Out of every side there 
'sprouts' two segments, giving 
rise to KS\. 



(d) G U KSi. The arrows in- 
dicate the cells Ci t k, 1 < k < 
3-4° = 3, of KSl ' 



Figure 19. An illustration of Definitions 14.111 and 14.121 in terms 
of KSq — A and KS\. The ghost of the segment So,k> denoted by 
go,k, is a middle-third segment of so.fc an d is removed from sq^ so 
that we may construct the cell of KS\. Then <?o.fc is referred 
to as the ghost of the cell C\,h of KS%. 

4.1. Compatible sequences of piecewise Fagnano orbits of fi(KS n ). In 

Q(KSo), the direction = 7r/3, measured relative to the fixed coordinate system 
described in Figure [T71 gives rise to two different orbits: The Fagnano orbit and 
a non- Fagnano orbit in the direction of 7r/3; see Figures I^TJl and The Fagnano 
orbit is the orbit in the direction 7r/3 starting at the midpoint of the base of A. 
This is, in fact, the shortest orbit of i7(A) (see jBaxUmj ). Any other orbit in the 
direction 7r/3 is necessarily twice as long as the Fagnano orbit & ' . 

We next discuss a generalization of the Fagnano orbit for the Koch snowflakc 
prefractal approximations il(KS n )- 

Definition 4.14 (Piecewise Fagnano orbits of fl(KS n )). If \p n € KS n is com- 
patible with a midpoint of a ghost g n -i,k of a cell C n ^ of fl(KS n ) in the 
direction 9(n/3) (which denotes an integer multiple of n/3 such that 0(n/3) is an 
inward pointing direction at y°), then the orbit ^„(x° , 9(n/3)) is called a piecewise 
Fagnano orbit of fl(KS n ). 

A piecewise Fagnano orbit ^ n (y°, #(7r/3)) is named as such for the fact that 
one can view such an orbit as the result of appending scale n copies of the Fagnano 
orbit of Sl(KSo) to every basepoint of an orbit ^ > „_i(y" 1 ,0(n/3)), where y^-i is 
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A 



Figure 20. The first orbit is the Fagnano orbit of the equilateral 
triangle billiard fl(KSa). The second orbit is the primary piecewise 
Fagnano orbit of Q(KSi) and the third orbit the primary piece- 
wise Fagnano orbit of fl(KS2)- These piecewise Fagnano orbits 
are called primary, because the first clement in the compatible se- 
quence of orbits to which all of these orbits belong is the Fagnano 
orbit of il(KSo); see also Definition 14. 161 

the midpoint of a ghost g n —i t k of a side s n -i t k of Q(KS n ) referred to in Definition 
14.141 see Figures l20l and I2T1 We denote a piecewise Fagnano orbit of Q(KS n ) by 

p&rC 1 , where y° is the initial basepoint of the orbit and 9(tt/3) is the initial inward 
pointing vector. When is collinear with x$ = 1/2 G I in the direction 9(ir/3) (or 
when, for some integer k n , a basepoint y^ n of the orbit & n {Vni 9{tt/3)) is collinear 
with Xq — 1/2 € / in the direction 9'(tt/3)), then we write the orbit & n {Vn , 9(n/3)) 
as pp^ n for the primary piecewise Fagnano orbit of fl(KS n ); see Figure [20l 

Remark 4.15. Later, we will see that for every orbit 9(ir/3)), there exists 

a unique element XqGI (the unit interval [0, 1] viewed as the base of A = KSq) 
such that Xq is compatible with a basepoint y%™ of the orbit n ^V%, , 0(t/3)) in the 
direction Att/3. (See Proposition 14.201 and the discussion preceding it.) Conse- 
quently, we will eventually write ^„(y", 9(w/3)) as i^„(a;° , tt/3), since i^„(a;" , rr/3) 
determines, and is determined by, the orbit ff n {y^, 9(n/3)). Then, in the case when 

^„(x°,7r/3) is a piecewise Fagnano orbit, we denote the orbit by p.^n° , since Xq 
determines the point and vice- versa. 

Definition 4.16 (Primary piecewise Fagnano orbit of Q(KS n )). Let p^n a be 

a piecewise Fagnano orbit of ft(KS n ). If Xq = 1/2 6 /, then p,^n° is called the 
primary piecewise Fagnano orbit of Q(KS n ). We denote the primary piecewise 
Fagnano of Q(KS n ) by pp^ n . 

We define a compatible sequence of piecewise Fagnano orbits as follows. 

Definition 4.17 (Compatible sequence of piecewise Fagnano orbits). The se- 
quence of orbits {&i{yi, 9(<fi))}°^ is called a compatible sequence of piecewise Fag- 
nano orbits if {&i(yf, 9((f)))}^ =0 * s a compatible sequence of periodic orbits (in the 
sense of Definition I4.8[) and there exists N > such that for every n > N, 
midpoint of a side s ni k of a cell C n ,k' of Q(KS n ) (1 < k' < 3-4 n ~ 1 ). Alternately, we 
say that a compatible sequence of periodic orbits is a compatible sequence of piece- 
wise Fagnano orbits if there exists N > 1 such that for every n > N, ff n (y^, #(</>)) 
is a piecewise Fagnano orbit of rt(KS n ); see Figure l20l 




28 



MICHEL L. LAPIDUS AND ROBERT G. NIEMEYER 



Figure 21. Examples of orbits in the direction ir/3. The first or- 
bit is a non Fagnano orbit of the equilateral triangle. The second 
orbit is a non piecewise Fagnano orbit with an initial direction of 
7r/3 (specifically, the second orbit will be an element is an eventu- 
ally constant compatible sequence of periodic orbits; see £14.21) . The 
third orbit shown in Cl(KS2) is an orbit for which the next orbit 
in the corresponding compatible sequence (to which the two orbits 
would belong) is a piecewise Fagnano orbit of V,(KSs), since each 
basepoint of the orbit lies on a midpoint of some ghost of a cell 
C 3 , k of Q(KS 3 ). 




We note that it follows from Theorem 14.61 that, under the assumptions of Def- 
inition OH {(y°,0(0)°)}£L o is a compatible sequence of initial conditions, in the 
sense of Definition 14.51 Moreover, it follows from Definition 14.171 that for every 
i > 0, there exists ki £ N such that 0(</>)° = kin/3. 

As alluded to above in Remark 14.151 we want to be able to describe every 

piecewise Fagnano orbit p^n n of Q(KS n ) in such a way that p^rC 1 is actually an 
element of some compatible sequence of piecewise Fagnano orbits determined by a 
particular orbit 0o(xq,tt/3) of fl(KSo), with x® an element of the base of KSq. In 

other words, if pJ^n™ is an orbit of ft(KS n ), then we want to show that there exists 
x„ £ KS n such that 1) midpoint of a side s n .k, 2) x° is collinear with x® 

and 3) (jMtt/3)) £ n (x° n ,n/3). 

We begin establishing such a connection by considering the following three 
contractive similarity transformations <f>i, fa, ^3 : R — > K: 

<j>i(x) 

(4.1) <f> 2 (x) 

Consider $ : K — > K, a map defined on the space K of all nonempty compact 
subsets of K, as $(K) := Uj=i fiji^)- When K is equipped with the Hausdorff 
metric, K becomes a complete metric space. (See, e.g., |Ba|, IEd|, IFaj .) Since each 
4>j is a contraction mapping, it follows from Hutchinson's Theorem [Hut that $ is 
a contraction mapping on K and by the Contraction Mapping Principle, $ has a 
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Figure 22. In this figure, we demonstrate where exactly the 
points of M(3) would be located in 1) Q(KS ) and 2) Cl(KS 3 ). 

unique fixed point attractor in K. If I = [0, 1], we see that $(/) = U*=i 4>j{I) = I> 
meaning that the unique fixed point attractor of $ is the unit interval I. 

For each integer n > 0, we define M(n) := $"({1/2}), where denotes the 
nth iterate of $ and $° := Id; see Figure Then it can be easily checked that 



(4.2) M(<*f) := |J M(n) 

n=0 

is the collection of all elements of the unit interval I with ternary expansions ter- 
minating in l's. Note that each M(n) is finite and hence, that M( c (a) is countably 
infinite. Furthermore, M(^) is dense in the unit interval /. 

The following lemma will help justify a key step in the proof of the next result 
(Proposition |4.19[) . 

Lemma 4.18. Let N > 1. For every 1 < i < N, there exists j% & {17 2,3} and 
Pj t € {0, 1, 2} such that 



(4.3) <p jN o tt> jN _, o ■ ■ ■ o 4> h (1/2) = Y, + J^n- 

i=i 

PROOF. We know that $ 1 ({l/2}) = {1/6,1/2,5/6}. More to the point, each 
number in $ 1 ({l/2}) may be given by §f + 5^" for a suitable choice of ji € {1, 2, 3} 
and p n € {0, 1, 2}. That is, 



r ^ = §+2^ ifji- 1,^=0 

<fe(l/2)= I 3 = 3 + 33 ifii = 2, ftl =l 

Next, let us proceed by induction. Let N > 1. Suppose that for each i = 
1, N, there exist jj e {1,2, 3} and € {0, 1, 2} such that 
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N 

*(i/2) = E Pl 



Then 



3 jv-i+i 2-3^' 
»=i 



f+iOfcw - ° 0ii(l/ 2 )) = 4>3n+i (^Z 3/-4+1 + o - ^) 

= / v Pji 4- 1 ^ I i 

Vfe 3^-'+! 2 • 3^ 3 3 

= V ^ + 

3(jv+i)-i+i ^ 2 • 3 Ar+1 ' 

4 = 1 



where p jN+1 G {0,1,2}. 

Therefore, for every TV > 1, we have 



AT 

n(V2) = E ^ 



3./V-4+I 1 2 • 3 N ' 

as desired. □ 

Proposition 4.19. Let x E I. Then x G M(^) i/ and only if £/iere exists a 
nonnegative integer N such that for every 1 < i < AT, f/iere exist 1 j< G {1,2, 3} and 
Pii G {0j 1) 2} suc/i i/iai 

Pit 1 



(4.4) x = 



i=i 



3 iv-i+i '2-3 



Proof. Let x € / such that a; G M(^). Then x G Af(iV) for some N > 0, 
where AT may be taken as the least such integer N . Since M(N) = $ w ({l/2}), we 
have that x £ $-^({1/2}) and, by Lemma [4.18[ for every 1 < i < N, there exist 
ji G {1, 2, 3} and p u G {0, 1, 2} such that 

N 1 

(4-5) a; = o o ... o ^(1/2) - ^ ^. +x + j-^. 

4 = 1 

Conversely, suppose now that there exists N > 1 such that for each 1 < i < N, 
there exist j 4 G {1,2, 3} and pj { G {0, 1, 2} such that 

N 1 



giV-i+l 2 . 3^ ' 

4=1 



Then, we see that for N — 1, 



l 



31-4+1 2 . 3 i 3 2 _3 

2—1 
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for some ji 6 {1, 2, 3} and pj 1 € {0, 1, 2}. Furthermore, if N > 2, we have succes- 
sively: 




3 J 3 3 



PJL + J_U 
3 2-3/ 3 " r 

= </> ijv o ... o 0^(1/2). 
This concludes the proof of Proposition 14.191 



P2,J 

3 



PjN 

3 



The number 5Z i=1 



Pii 



□ 



3W _ i+1 is an endpoint of an interval A C I of length 1 /3 . 
By the very nature of the Koch snowflake construction (via an iterated function 
system), there is a side SN,k of the prefractal approximation D,(KSn) such that the 
interval A is a translate (or a translate and rotation by ±7r/3) of this side sn^ in 
the direction 4ir/3. 

Let Xq <E M(N). Then, by Proposition 14.191 there exists {pj t }fLn with ji € 
{1,2,3} andp l: , € {0,1,2} for each i = 1,...,N, such that xg = Y^Li 



3 JV-i + l 



1 



meaning, Xq is the midpoint of the interval in the base of the equilateral triangle for 



which 3 N P -i+i is an endpoint. Therefore, in the direction 47r/3, the endpoint 

of the line segment connecting Xq with the segment SN,k of the boundary KSn 
is a midpoint of the segment SN,k- Denote this point by x N . Since there are 
no points of Q(KSn) \ on the aforementioned line segment, we know that 

we may construct a compatible sequence of initial basepoints {x®}^ such that 
{&i(x®, 7r/3)}^ is a compatible sequence of piecewise Fagnano orbits. 

Conversely, we want to show that if is a midpoint of a side s n .k of Q(KS n ), 
then, for some N > 0, & n {Vn, 0(t/ 3 )) determines an element x% G M(N) C M(<*f), 
and ultimately, a compatible sequence of initial conditions that corresponds to a 
compatible sequence of piecewise Fagnano orbits {£?i(x®, 7r/3)}°^ |3 Consider an 
inward pointing direction 6(ir/3) at a midpoint y° £ if on a side s ni fc of ^(ifSVi)- 
The billiard ball must reflect off of at most one other side of a cell in which the side 
s n ,k was a part of (this may not be a cell C n ^, but rather a cell C n /^k', with n' < n 
and fc' < 3 • 4" ) before exiting the cell (by convention, a cell Co,fc is the equilateral 
triangle); see Figure [23] 

By means of the local and global symmetry of the prefractal billiard table 
Q(KS n ), as shown in Figure [231 and discussed in the corresponding caption. One 
can deduce that there must be a basepoint y^™ of the orbit ff n (y^, 9(ir/3)) that 
corresponds to a midpoint and that, after colliding with the boundary KS n , the 



^Recall that the angle 0(tt/3) is an integer multiple of 7r/3, measured with respect to the 
fixed coordinate system shown in Figure [TT] 
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Figure 23. Depicted in this figure is a classical tool for analyzing 
billiard orbits: unfolding the orbit. We consider a billiard orbit's 
trajectory entering a cell C„ at an angle of 7r/3. As expected, the 
billiard ball must exit the cell after two reflections and at an angle 
of — it/3. We verify this by considering the unfolded trajectory in 
the plane and noticing that the trajectory continues on unimpeded 
as it passes through a reflected copy of the opening. (See Figure |2"41 
for a generalization of this discussion to the case of two collisions 
with the boundary.) 

billiard ball next traverses the interior of £l(KS n ) H fl(KSo) in the direction of ir/3. 
We claim that upon doing so, the billiard ball must pass through some element of 
Mi^to). Indeed, if it did not, then could not have been a midpoint of any sideF^ 
Specifically, y° has a value g^pr, relative to s n ,k, and y^ n has same the same value, 
but relative to some side s n ,k', 1 < k' < 3 • 4™. Since s n ,k> can be translated (or 
translated and rotated by ±7r/3) in the direction of 7r/3 so as to correspond with 
an interval of length 1/3™ in /, it follows that the billiard ball passes through an 
element of M(^); call this element Xq. Then, if we define x^ by a;" := y^ n , the 
orbit G n (x Q n , 7r/3) contains (y^,#(7r/3)) and constitutes an element in a compatible 
sequence of piecewise Fagnano orbits. 

We summarize our discussion in the following proposition: 



■^Recall that y„ n can also be viewed as the image of y„ under the action of particular local 
and global symmetries of Q(KS n )- 
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Figure 24. This figure generalizes what we have seen in Figure 
l23l An orbit's trajectory that enters a cell C n parallel to a side of 
C n will exit after two collisions. 

Proposition 4.20. Let n > and 0(n/3)) be compatible with a piecewise 

Fagnano orbit p<^jy V • Then there exists a unique element Xq G M^) and a com- 
patible sequence of piecewise Fagnano orbits {ff^xf, 7r/3)}££ such that 

(4-6) & n (y n ,6(<ir/3)) € {^(x°^/3)}£ 

and 

(4.7) p&fr = ff m {x Q ml tt/3), /or every m > N. 

Conversely, for every Xq € M(&), there exists n > and a piecewise Fag- 
nano orbit p^n n of Q(KS n ) such that x^ is compatible with Xq and &q(xq, 7t/3) 

x° 

is compatible with p^ n ™ • Moreover, there is a least integer n > such that x n is 
compatible with Xq and for every integer k > n, the compatible orbit €?k(x1,9®) is 

a piecewise Fagnano orbit p^ n " determined by Xq . 

We then are in a position to say that Xq £ Mifio) determines a piecewise 

x° x° 

Fa gnano orbit p^ n n and, conversely, pJ^n™ determines a unicjue element Xq £ 

x° x° 

M\€), allowing us to write p^ n n as p^ n ° without any ambiguity. 



34 



MICHEL L. LAPIDUS AND ROBERT G. NIEMEYER 




Figure 25. For n > 1, the symmetry group of KS n is the dihedral 
group D 6 (see fc|2.2[) . Local symmetry is seen at the level of a cell 
C n of Q(KS n ). The symmetry of interest is symmetry with respect 
to the angle bisector of the acute angle of the cell C n . 

4.2. Eventually constant compatible sequences of periodic orbits. 

Consider the ternary Cantor set ff. It is the collection of all elements in the unit 
interval I such that an element has a ternary expansion consisting only of O's and 
2's. This is not to imply that an element of c € will not have an expansion consisting 
of any number of l's, but that such a number must have some ternary representa- 
tion satisfying this rule. For example, denoting a base-3 number in I by O.M1U2...3, 
we see that 1/3 = O.I3 is in ^ because 1/3 = 0.023, which has no l's. Similarly, 
1/4 = O.O23 is also in . However, 1/2 = O.I3 is not in because 1/2 does not 
have an equivalent representation that would satisfy the rule for being in 

^E3 The 

complement c rf c — I \ ^ (with respect to the unit interval I) can be described as 
the collection of elements in I having a ternary expansion comprised of at least two 
l's. For example, 0.1 13 is not in the Cantor set, because it lies in the open interval 
(1/3, 2/3), which is removed during the geometric construction process as shown in 
Figure H (See Example [2T0] of 32~3l l 

Let be the collection of all elements of ^ with infinite ternary expansions 
that do not terminate solely in a single character (repeated ad infinitum). These 



21 Indeed, 1/2 is in the deleted interval (1/3, 2/3) at the first stage of the construction of V. 



FAMILIES OF PERIODIC ORBITS OF THE KOCH SNOWFLAKE FRACTAL BILLIARD 35 

elements are naturally represented as infinite sequences of O's and 2's, with peri- 
odic expansions representing rational values and aperiodic expansions representing 
irrational values of ^ . 

In any prefractal approximation ^ n of every x £ < £' is contained in a (one- 
dimensional) connected open neighborhood U x such that UxD^n = U x . Therefore, 
for every n > and x® g c €' C 7, ^„(x" , 7r/3) is a well-defined periodic orbit of 
fi(i ; £f<S'„)I3 Moreover, the compatible sequence of periodic orbits {^(2;°, 7r/3)}°^ 
is an eventually constant sequence of orbits, since x® £ KS n , for every n > 0. 

We may extend this result. In order to do so, we first recall Proposition 
14.201 Xq € Mffi) determines and is determined by a piecewise Fagnano orbit 
of Q(KS n ), n > 1. The proof of this statement (given in the discussion preceding 
the proposition) relies heavily on the local and global symmetry of the prefractal 
snowflake tt(KS n ). We take advantage of such symmetries when constructing an 
orbit G n (xP n , 7r/3) of ft(KS n ) with an initial basepoint on the boundary KS n of 
n(KS n ). 

Consider ^„ : fc C S n ,k a scaled copy of where n > and 1 < fc < 3 • 4™. 
Then, define fc to be the collection of all points of c £ n ,k with a ternary expansion 
consisting of infinitely many O's and 2's. If y° £ ^n.k' then the orbit G n (^, 8(n/3)) 
is a well-defined orbit of Q(K S n ). In addition, there exists a basepoint y^ n of the 
orbit ff n {y Q n , 9(n/3)) such that after colliding in the boundary at y^ 1 , the billiard 
ball then traverses the interior of £l(K S n ) nfl(KSo). In analogy with the case when 
y^ n was a midpoint of a side s n ,k, the billiard ball must now pass through a point 
Xq € I with a ternary expansion consisting of finitely many l's and infinitely many 
O's and 2's. Moreover, this compatible sequence is eventually constant, meaning 
that there is a positive integer N such that G n {xP n , 7r/3) = &^(xP N , 7r/3), for all 
n > N. We have therefore proved the following proposition. 

Proposition 4.21. Let N > and k < 3 • 4 . If y% £ k' then there exists 
x° N £ KSn compatible with some element Xq £ I having a ternary representa- 
tion consisting of finitely many 1 's and infinitely many 's and 2 's and such that 
{0i(x®, 7r/3)}[f is an eventually constant compatible sequence of periodic orbits with 

Remark 4.22. If there exists a positive integer N > Osuchthat {@i{x\, 7r/3)}i>Ar 
is a constant compatible sequence of periodic orbits, then {^(2;°, 7r/3)}°^ is an 
eventually constant compatible sequence of periodic orbits. 

As a result of Proposition 14. 2 1 1 and Remark 14.221 we have the following defini- 
tion. 

Definition 4.23. If Xq £ I has a ternary expansion consisting of finitely many 
l's and infinitely many O's and 2's, then we call the resulting compatible sequence 
of periodic orbits {&i(x\, 7r/3)}°^ an eventually constant compatible sequence of 
periodic orbits. 

Example 4.24. Consider the element Xq £ I with the ternary representation given 
by O.IO2O3, where, as usual, the overbar indicates that the corresponding string 20 is 
repeated ad infinitum. This is the value 1/3 + 1/12 = 5/12, which is not an element 
of Mffi). Then, choose x® in KS\ on a side s\ y k such that x\ is collinear with Xq 



'This means that it does not hit any of the vertices of KS, 
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in the direction 7r/3, relative to the fixed coordinate system. This point identified is 
the point O.O23 scaled by 3, or 1/4 scaled by 3, residing on the side si^. Since 1/4 
is a point of the Cantor set (because 1/4 = 0.02 3 ), 1/12 is a point of the Cantor set, 
and hence, the point x\ is the point 1/12 on the side and in the scaled Cantor 
set^ n . k B 

Consequently, the orbit ofQ(KSi) is compatible with that offl(KSo) and the 
corresponding compatible sequence of periodic orbits is an eventually constant com- 
patible sequence of periodic orbits. 

If Xq is an element with a ternary expansion consisting of infinitely many O's 
and 2's and finitely many l's, then, by Definition 14.231 the compatible sequence of 
periodic orbits {0i(x®, 7r/3)}°^ is an eventually constant compatible sequence of 
periodic orbits. If N > is such that {^(x® ,ir/3)}°l N is a constant sequence of 
compatible orbits, then, for each n > TV, we call the corresponding orbit G n {x a n , tt/3) 
a ^ -orbit of il(KS n ). In the context of the fractal billiard Sl(KS), we will refer to 
a ^-orbit as a stabilizing periodic orbit (or, simply, a stabilizing orbit) of the Koch 
snowflake billiard Q(KS); see M5.41 

4.3. Compatible sequences of generalized piecewise Fagnano orbits 
of J7(KS n ). A recurring theme thus far in §4 is that if 8(ir/3)) is an ini- 
tial condition of a periodic orbit ^„(y", 9(n/3)) in Sl(KS n ), then there exists 
[(a£,7r/3)] e (niKS^xS 1 )/ - such that (x°,tt/3) e ^„(y°, 0(tt/3)), x° n iscollinear 
with a point x[J in the unit interval I and {<^(x°, 7r/3)}^ forms a compatible se- 
quence of periodic orbits. The points x of the billiard table Q(KS n ) for which we 
have shown this to be true are 

(1) x £ fl(KS n ) such that x <E <£' k C s n . kl for some k < 3 • 4"; 

(2) x € J^lfSn) such that for every n' > n, there exist fc' < 3 • 4" and 
x' € ^/ fe/ C Sn'./c' such that x and x' are compatible in the direction of 
0(tt/3); ' 

(3) x € s nj fe such that x € M(f£) n ^, the set M(f£) n ^k being a scaled (or 
scaled and rotated) copy of M(^), viewed as a subset of s n ^- 

A point y° on a side s ni fc not falling into any of these categories (and which, for 
any m > 0, is not compatible with a corner of ft(KS m ) in the direction of 9(tt/3)) 
constitutes an initial basepoint for which we call the resulting orbit & n {Vni 
a generalized piecewise Fagnano orbit of £l(KS n ). Because we are measuring angles 
with respect to the fixed coordinate system, recall that 9(ir/3) is our convention 
for indicating that 7r/3 was the measure of the angle when measured relative to 
the side s n> fc on which y° lies and 6(jr/3) is then the measure of the same angle 
measured relative to our fixed coordinate system. We note that such an orbit is 
never a piecewise Fagnano orbit nor is it an orbit ^/(y^,, 6{tt/3)) such that y^ 
and are collincar in the direction 8(ir/3) and £? n /(y°/, &(tt/3)) is a piecewise 
Fagnano orbit for any n! > n. As was done before, we may find x° € s n ,fc such 
that <^„(x° , 7r/3) = ff n {y^, 6(ir/3)) and x° is collinear with some x[j G J so that 
{^i(x", tt/3)}°^ constitutes a compatible sequence of periodic orbits. We call 
such a compatible sequence a compatible sequence of generalized piecewise Fagnano 
orbits. 



''This can be more clearly seen if one views the side Si j as a rotation and translation of the 
left middle third of the unit interval / (viewed as the base of KSo). 
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The element Xg in the corresponding compatible sequence of basepoints is not 
an element of M( c (o)\ otherwise, there would exist N > and & n {x% , it / 3) in 
the compatible sequence such that the basepoints of &n (x%, 7r/3) would have to 
be midpoints of particular sides of some prefractal approximation Q(KSisr). In 
addition, if x® N is the initial basepoint of a generalized piecewise Fagnano orbit, 
then the initial basepoint Xg £ I of a compatible orbit &q(xq, 7t/3) of Q,(KSq) has a 
ternary expansion consisting of infinitely many O's and l's, 2's and l's or O's, l's and 
2's; otherwise, ff^/^x^ , 3 ,) would be an orbit in an eventually constant compatible 
sequence of orbits or a compatible sequence of piecewise Fagnano orbits. 

4.4. Properties of periodic orbits of f2(KS n ) in the direction of tt/3. 
According to the preceding discussion, in each prefractal approximations Q,(KS n ), 
we may collect orbits according to the nature of the basepoints, or as piecewise Fag- 
nano orbits, %f '-orbits and generalized piecewise Fagnano orbits. Note that this does 
not determine equivalence classes of orbits in the direction tt/3. Indeed, if x° and y° 
are two elements of a side s Uy k of Q(KS n ), then n {x% , 0(ir/3)) and @n(ifl , 6(ir/3)) 
will be equivalent orbits, i.e., they will have exactly the same lengtho It is only 
when x° and y° differ in their nature, e.g., x° £ ^ k and y° £ ff^k, that we see a 
difference in the orbits in each compatible sequence given to us by [fi^x®, tt/3)}°Z q 
and {@i{yi, 7t/3)}°^ - The significance of this fact will become apparent in §5 where 
it will aid us in giving a description of periodic orbits in the direction of 7r/3 of the 
Koch snowflake billiard Q(KS). 

Remark 4.25. Let n > 1. If & n {x® n , 7r/3) and <^„(y°, tt/3) are two periodic 
orbits of Q(KS n ) with the same period (that is, #<^„(x° , tt/3) = #<^n(y° , tt/3)) 
such that, for every k n < #<^„(x° , #°), x^" and y^™ lie on the same side of Q(K S n ), 
then the corresponding paths traced out by connecting consecutive basepoints have 
exactly the same length. Such a fact follows from the known equivalence between 
the billiard flow on the rational billiard ft(KS n ) and the geodesic flow on the 
associated prefractal flat surface S(KS n ); see the discussion at the very end of 

Let n > 1. If (x[j)j denotes the ith character in the ternary expansion of Xg (in 
terms of 0, 1, 2), then define u;„(xq) to be the cardinality of the set {(xq)^ | (xg) 4 = 
1, 1 < i < n}. 

Proposition 4.26. Let n > 0, fc < 3 • 4™ and x° <E s n ^- Then there exist m < n, 
h < 3 ■ 4 m and y^ £ s m ^h such that ff m (y^ n ,Tr/3) is a piecewise Fagnano orbit of 
fl(KS m ) and y^ compatible with some y® £ M(^) with the following being true: 
there exists (y^ 1 ,0(w/3) £ ^' m (y^ n ,n/3) such that x^ and y 1 ^ 1 both lie on the same 
side Sm^h' of Q,{KS m ), h' < 3 ■ 4 m ; and the ternary representation o/x° and y^ 
(relative to s m .h>) are identical up to the first m — 1 characters. 

PROOF. We defer the proof of this statement until the very end of £)5.2.U since 
much of the machinery introduced in §5 makes proving this proposition considerably 
less tedious. □ 

Theorem 4.27 (Computation of the period of n (x° n , tt/3)). Let {^(x?, tt/3)}£ q 
be a compatible sequence of periodic orbits. Then every orbit G n (xP n , 7r/3) in the 
compatible sequence has a period #&n{Xn, 7r /3) that is determined by the number 



See i|2.1l for a more thorough discussion of the equivalence relation on orbits. 
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(d n (xo)- Specifically, for every n > 0, the orbit #^„(x,j, 7r/3) is given by the fol- 
lowing formula: 

(4.8) #<?„(s°,7r/3)=3-2 w »<*8>. 

PROOF. Let {@i{x\, 7r/3)}°^ be a compatible sequence of periodic orbits. For 
the purpose of readibility, we consider the following three cases: 

(1) {i^i(x?,7r/3)}°^ is a compatible sequence of piecewise Fagnano orbits. 

(2) {0i(x®,n/3)}°Z o is an eventually constant compatible sequence of periodic 
orbits. 

(3) {&i(x®, 7t/3)}°^ is a compatible sequence of generalized piecewise Fag- 
nano orbits. 

Case 1: Suppose {£?i(x°, tt/3)}^ q is a compatible sequence of piecewise Fagnano 
orbits. Then there exists N > such that for every n> N, & n {xP n , tt/3) is a piece- 
wise Fagnano orbit and for every n < N, @ n {x® n , 7r/3) is not a piecewise Fagnano 
orbit. Let n = 0. Then ^ n (a;° , tt/3) has period #^„(x° , 7r/3) = 3 if the orbit is 
the Fagnano orbit of S}(KSq), and #^„(a;° , 7r/3) = 6 otherwise. 

Now suppose that n = 1. Without loss of generality, we may assume 2q <E 
(1/3,2/3). Therefore, Wi(xg) = 1. In addition, we see that $^ n (x^, 7r/3) = 6 = 
3 . 2 Wl ( a: o). 

Next, let M > 0. Suppose that for every n < M, we have #^ n (a;° , tt/3) = 
3 . 2 w "( a: o). We may as well assume that N < M. Therefore, e M {x° M ,^/^) is a 
piecewise Fagnano orbit of H(KSm)- By definition, p&m+i is constructed from 
P-^m appending 3 • 2 WM ^ X ° S > many scaled copies of ppJ^o (the Fagnano orbit of 
the equilateral triangle billiard Sl(KSo)) to each basepoint of p^m ■ Hence, 

(4.9) = 2 ■ #p&$ = 2 • 3 • 2 Wm(x « ) = 3 ■ 2 UM ^ +1 . 

Since (xq)„ = 1 for every n > N, we deduce that wm+i^o) = wm(^o) + 1> an d the 
result follows for a compatible sequence of piecewise Fagnano orbits. 

Case 2: Suppose 7r/3)}^. is an eventually constant compatible sequence 

of periodic orbits. Then there exists N > such that for every n > N, (? n , tt/3) 
is identical to &n{x° n , tt/3). Moreover, u n (i5) = ojjv^o) f° r every n> N. Suppose 
iV = 0. Then #0 N (x%,n/3) = 6. Suppose N > 0. By PropositionH^H there exist 

M < N and y® E I such that pJ^T? forms a piecewise Fagnano orbit of VL{KSm), 
u>n(xq) = oJm(Vo) an d #^ > A'( a; Ar: t/3) = #P=^f . In Case 1, we saw that the period 
of a piecewise Fagnano orbit was determined by the formula #p&$ = 3 ■ 2 UM< ^°) , 
Therefore, 

(4.10) # € N (x%, tt/3) = #p&$ = 3 ■ 2 WM ^) = 3 ■ 2 6W <" ! °>. 

Since G n (xP n , tt/3) is identical to /?n(x%,tt/3), the result follows for an eventually 
compatible sequence of periodic orbits. 
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Case 3: Suppose {&i(x\, 7r/3)}££ is a compatible sequence of generalized piecewise 
Fagnano orbits. For each n > 0, there exists Vq £ I and M < n such that v^ V m ^ s 
a piecewise Fagnano orbit of SI(KSm) and w n (a;g) = u)m{Vq)- 

Let N > and suppose that for every n < N, #ff n (x° n ,ir/3) = 3 • 2^^o\ 
Then, there exist v.® £ I and M < N + 1 such that v&'m ^ s a piecewise Fagnano 
orbit of VI(KSm) and u>jv+i(2;g) = ujm(Vq). Therefore, 

(4.11) #^,v +1 (x? v+1 ,7r/3) = = 3 • 2 w «^o°) = 3 . 2 "w+i<»8). 

This concludes the proof of Theorem 14.271 □ 

Notation 4.28. We now explain the notation that is about to be used in the fol- 
lowing theorem and proof. The characteristic function x *s defined on the space of 
characters {0, 1,2} and is given by 



(4.12) x [a] :-- 



if a = 0, 2 

1 if a = 1 ; 

in other words, it is the characteristic function of {1}. A ternary expansion of an 
element Xq £ I begins with 0. • ■ ■ . Consequently, in order to simplify the notation, 
when we discuss x[( x o)n]i we are viewing the ternary expansion of Xq as the sequence 
of characters occurring to the right of the decimal point (and we no longer indicate 
the subscript 3 in the expansion) . 

Theorem 4.29 (Length of the billiard ball path corresponding to ff n (x^ , tt/3)). 
U {&i( x i 7 7r /3)}^o * s a compatible sequence of periodic orbits and Jzf is the length 
of the Fagnano orbit of the equilateral triangle billiard il(KSo) (i.e., the shortest 
inscribed polygon in ^I(KSq)), then, for each n > 0, the length \@ n (x® n , 7r/3)| of 
the path traced out by connecting consecutive basepoints of the orbit & n (x Q n , 7r/3) is 
given by 



(4.13) |^ n (4,V3)|=2if + ^x[(^)J#^-i(^i^/3) ¥ , 

i=2 

where, as before, (xq), denotes the ith character in the ternary expansion represent- 
ing xg. 

Proof. Let n = 0. If 1 is the ternary expansion (in terms of the characters 
0, 1, 2) of an element Xq £ I, where I is viewed as the base of KSo, then &q(xq, 
is the Fagnano orbit of H(KSq) and has length | (xg , 7r/3) | — ^ ■ If nas 
a ternary representation different from 1, but beginning in the character 1, then 
|^b( a; o> 7r /3)| = 2 Jzf. In either case, if x\ is collinear with xg, then \€?i(x1,ir/3)\ = 

Consider the basic case n = 2. Let x^ £ KS2 be collinear with Xq £ I (viewed 
as the base of the equilateral triangle KSq). Then, x\ is the basepoint in the com- 
patible orbit &2(x%,k /i). Let (xgL be the ith. character in the ternary expansion 
of Xq. We want to show that 

2 

(4.14) 1^(^,^/3)1 =2^ + ^ X [(^)J#^-i(^_i, 7 r/3)f. 
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If (a;g) 2 = 1, then there exists yg £ M{<#) such that \0 2 (x%,n/3)\ = \p^ \ and 
(xo) i = (j/o)i f° r every i < n = 2. There exists yj 1 G KS\ that is collinear with 
yg such that \p&%°\ = \&i(y%, tt/3)| + (y?, tt/3)^, since p^f is a piecewise 
Fagnano orbit of KS 2 - Since yg and xg have the same first two characters in their 
respective ternary expansions, we have that 



1^(^^/3)1 = 1^(^,77/3)1 

and 

#^ 1 (j/ 1 , 7 r/3) = #^ 1 (^,7r/3). 

Since <^i(yf,7r/3) = 2Jzf (which is independent of the choice of basepoint, so 
long as such a choice is not a corner of the billiard table Sl(KSi)), it follows that 



(4.15) \0 2 (xl tt/3)| = 2j2f + £ #^_ 1 (x°_ 1 ,7r/3)^. 



i=2 



If (zg) 2 ^ 1, then X [(zg) 2 ] = and 



(4.16) |^ 2 (z°, tt/3)| = |^(x?,7r/3)| = 2JS?. 

In either case, we have shown that Equation (|4.14|) holds. 

Let us now proceed by induction and fix TV > 2. Suppose that for every n < N , 



(4.17) |^„(4,7r/3)|=2^ + ^ X [(x°) i ]#^_ 1 ( a; ti,7r/3)- r . 

Then there exists yg G M(^) and M < iV + 1 such that for all i < M, (a:g) i = (yg) € 
and 



(4.18) |^ +1 (x5V +1 ,7r/3)| = |p^|. 



If M = N + 1, then the nature of dictates that 



(4-19) Wt\ x \ = |^(^, 7 r/3)| + #^(y° r ,7r/3)^ T . 

Applying this fact, the induction hypothesis (|4.17j) and Remark l4.25l we have that 
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\0 N+1 (x° N+1 ,ir/3)\ - \pK°+i\ 

5? 



= \e N (y o N ,n/3)\+#0 N (y o N ,Tr/3) 
= \& N (x° N ,ir/3)\+#0 N (x° N ,Tr/3) 



3 jv+i 
'3— 

N 



(4.20) = 2^ + J2x[(x° ) i ]^ i -i(xl 1 ,n/3)§ 



i=2 



#^v(x° v ,7r/3)- 3N . , 

2V+1 

2Jgf+£ X [(a:8)J#^-i(a:?-i,T/3) 



i=2 



If M < TV + 1, then, again applying the previously mentioned fact, the induction 
hypothesis and Remark 14.251 we deduce that 



|^ +1 « +1 ,tt/3)| = \p^l°\ 

= |^f-i(ySf_i,7r/3)|+#/J' Jltf -i(tfS f _ 1 ,7r/3)^ 

= | G M - 1 (4,-1 , 7T/3) I + #<?M- ! , ^ 

M-l ~> 

(4-21) = 2JSf+X)x[(A]#^-i(^S-i^/3)3T + 

#^M-i«_ l5 7r/3)^ 
= 2j2f + 2x[(^) i ]#^-i(^-i,7r/3)- r 

z=2 6 
N+l „ 

= 2^+^ X [(x°)J#^_ 1 (^_ 1 ,V3)-, 

t =2 15 

where the last lines of the calculation in Equation (|4.21[) follow from the fact that 
the characters (x§) -, with M < i < iV + 1, are necessarily never equal to 1, meaning 
that x(x°)i = for M < i < N + 1. 

□ 



5. Periodic orbits of fi(KS) in the direction of 7r/3 

In §4, we were able to group orbits with initial directions 0(tt/3) into par- 
ticular categories: piecewise Fagnano orbits, "^-orbits and generalized piecewise 
Fagnano orbits. We stressed that this grouping was not equivalent to the classi- 
fication of orbits described in £12-11 Such a grouping is, however, meant to allude 
to a description of orbits in tt(KS) with an initial direction of 7r/3, where this de- 
scription is determined by the ternary expansion of an initial basepoint of an orbit 
0q(xq, it/ 3) with an initial direction of tt/3. In addition, we showed that for every 
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orbit ffniVni ^( 7r /3))) there existed x° £ KS n and Xg £ I collinear in the direction 
of tt/3 such that ^„(y°, 6 (tt/3)) = <? n (x°, tt/3). 

We want to somehow overcome the fact that the boundary of tt(KS) is nondif- 
ferentiable. When determining orbits of the snowflake, we want to stress the fact 
that we are primarily interested in the collision points with the boundary KS. In 
the theory of mathematical billiards, one may consider the billiard orbit to be the 
path traced out by the billiard ball or just the (equivalence classes of) ordered pairs 
(x, ff) £ (B x S 1 )/ ~. Even more simply, one may think of the path traced out by 
the billiard ball or the basepoints on the boundary B as a set of dynamically ordered 
points (3 When we want to make the distinction between the path and the ordered 
pairs in the prefractal billiard tables Q(KS n ) and in Ci(KS), we will explicitly refer 
to such a path as the billiard ball path (or, simply, the orbit) corresponding to a 
particular orbit and to the collection of collision points as the footprint of the orbit. 
In other words, the footprint is the intersection of the orbit with the boundary of 
the billiard table. Strictly speaking, we cannot speak of basepoints here, because 
this language is indicative of the existence of a well-defined phase space (a unit 
tangent bundle), which we have yet to rigorously establish. 

As noted in JQ in the sequel, we make a slight abuse of language. A self-similar 
set is the unique fixed point attractor of a particular iterated function system; see, 
for example, Figure [2] of fjl] By abuse of language, we will also say that a set 
is self-similar if it is the union of finitely many isometric, abutting copies of a 
given self-similar set, much as the snowflake curve is the union of three (isometric, 
abutting) copies of the Koch curve; see Figure [3] of £Q] 

We next define a "self-similar orbit." 

Definition 5.1 (Self-similar orbit). Let 7 be a periodic orbit of Q,(KS). Then, 7 
is said to be a self-similar orbit if its footprint is a self-similar subset of KS C M 2 . 

We want to show that it may be possible to appropriately recast particular 
compatible sequences of periodic orbits in such a way that a suitable notion of 
limit may be applied so as to yield an orbit of the Koch snowflake billiard ft(KS). 
Specifically, we will provide a plausibility argument as to why a suitable limit of a 
compatible sequence of piecewise Fagnano orbits constitutes an orbit of the Koch 
snowflake billiard tt(KS). In addition to this, we will show how an eventually 
constant compatible sequence of periodic orbits can be thought of as a periodic 
orbit of ft(KS), called a stabilizing orbit (or a ^-orbit) of Sl(KS). We will also 
discuss a connection between the ternary representation of the initial basepoint of 
the initial orbit of the eventually constant compatible sequence of orbits and the 
prefractal approximation fl(KS n ) at which such a compatible sequence stabilizes. 

5.1. Addressing system for KS. We now describe an addressing system for 
the Koch snowflake KS that will be used in the sequel. This addressing system was 
introduced (and used) by M. Pang and the first author in jLaPaj in their study of 
the boundary behavior of the eigenfunctions of the Koch snowflake drum. 

First, we mention that the ternary representation of an element of some side 
s n: k of Q(KS n ) will be given in terms of the alphabet {Z,c, r} (for left, center and 
right) so as to not confuse the reader later when we discuss the addresses of points 



In the theory of dynamical systems, such a set is called the Poincare section of the flow; 
see §2.2. 
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of the Koch snowflake curve KS in terms of the alphabet {0, 1, 2, 3, 4, 5}@ In the 
Koch snowflake curve, there are three types of points, as described below: 

(1) Corners: a vertex of the polygon KS n is a corner. As expected, such a 
corner in KS n is also a corner in KS m for all m > n. 

(2) Non-ternary points: an clement with a ternary representation (relative to 
a side s„,fe) not terminating solely in I's or r's. As such, a non-ternary 
point of KS is an element of some ^ k C s n ,fe, where ^ fc was defined in 
the paragraph preceding Proposition 14.211 

(3) Elusive limit points: an element of KS that is never an element of a 
finite approximation KS n , for any n > 0. Though every point of KS is a 
limit point, we make this abuse of language so as to emphasize the elusive 
nature of these points which are not vertices of KS n , for any n, and are 
not elements of c €' n k , for any n and any k < 3 • 4™ . 

In order to understand how these points are assigned addresses, we need to 
understand exactly how the addressing system works. The alphabet from which we 
will draw upon to construct addresses of points of the snowflake is {0, 1,2,3, 4, 5}. 
We start by noting that every side of a prefractal billiard table may be identified 
by a finite address. An address of a side may end in any character 0, 1, 2, 3, or 4 
but never 5. However, an address may start with 5. If an address of n characters 
ends in a 1 or a 3, then the side s n ,k identified by this address is not a subset of any 
of the sides of KS n -\. If an address of n characters ends in 0, 2 or 4, then the side 
s n> k is either the left third or right third interval of a side s n -i,fc' of S n —i), 
and has a length of 1/3™. In Figure [26l we label all the sides of £l(KSi) according 
to this convention. 

Heuristically, the labeling is such that the next character in an address is de- 
cided by 'moving forward' or shifting left or right. When moving forward, there is 
an option of a 'twisting' right or 'twisting' left, which results in appending a 1 or 3 
to the address, respectively. When making a shift to the left third or the right third 
of a particular side and the last character of an address is a 1, then the address 
is appended with a 4 or 2, respectively. Likewise, having an address end in 3 and 
performing a shift to the left third or right third of a side results in appending this 
address with a 2 or 0, respectively. Note that when 'moving forward' into a cell 
emanating from the side so,k with the label 5, one must first rotate the directional 
by 7r so that the proper addresses may be applied to the sides of the cell C n> fc; see 
Figure Ell 

In terms of the addressing system described above and in Figures [26j [27] and 
[28} the three types of points of KS listed above have very specific symbolic repre- 
sentations: 

(1) Points of the snowflake that are corners (i.e., vertices of some prefractal 
approximation KS n , with n > 0) have finite (i.e., terminating) addresses; 
cquivalently, the addresses of these points may be infinite, terminating 
solely in 0's, 2's or 4's. 

(2) Points that are elements of a set k when viewed as a subset of some 
side s n ^k of fl(KS n ) have addresses ending solely in cither 0's and 2's or 
solely in 4's and 2's. 



In §4, we used the standard alphabet {0, 1,2} to describe the ternary representation of an 
element in I. 
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Figure 26. The sides of the prefractal billiard il(KS{) labeled 
according to the particular addressing system described in M5. 11 




Figure 27. The directional used to identify points on the Koch 
snowflake K S. 

(3) Elusive limit points of the Koch snowflake KS have addresses given by 
(finitely or infinitely many) O's or 2's and/or 4's or 2's interspersed through- 
out a word consisting of: 

(a) infinitely many l's and 3's, 

(b) infinitely many l's and finitely many 3's, 

(c) or infinitely many 3's and finitely many l's. 

Example 5.2. The address 13131000220 identifies an element of ^ k on the side 
13131 of the prefractal Q(K S±) (and hence, a type (2) point). Moreover, the address 
131112343 identifies an elusive limit point of the snowflake curve KS (and hence, 
a type (3) point). Finally, the address 54 identifies the point 1/3 on the base of 
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Figure 28. Addresses of points of the Koch snowflake are not nec- 
essarily unique. In the context of a finite approximation KS n , the 
addressing system serves to label sides of the prefractal approxima- 
tion KS n - Points A, B and D then lie on sides which are labeled 
by the addresses 134, 131, 511, respectively. The point C, however, 
does not lie exactly on a side, but at a vertex. Hence, one could 
say that C lies on the side with address 504 or 534. Moreover, in 
the limit, C may be represented by the equivalent addresses 504 
and 534, where (as before) the overbar indicates that 4 is repeated 
ad infinitum. Since a billiard orbit is a collection of dynamically 
ordered points, it follows that each basepoint of an element of an 
orbit can be given an address in terms of this addressing system 
and still be distinguished from other basepoints of elements of the 
orbit with the same address. 

the equilateral triangle billiard Q(KSo); this is a type (1) point which will end up 
being an obtuse corner in every subsequent prefractal billiard table. 

5.2. The symbolic representation of a compatible sequence of initial 
basepoints. In 21 we defined what it meant to be a compatible sequence of orbits 
and a compatible sequence of initial conditions (see Definition l4.4l and Definition ^. 51 
respectively, along with Theorem l4.6p . Such a construction amounted to picking an 
initial basepoint Xq of an orbit &o(xq, 6®) such that the initial basepoint of every 
subsequent orbit in the compatible sequence was collinear with Xq in the same 
initial direction relative to a fixed coordinate system. 

Now that we have a clearly defined addressing system for points of the Koch 
snowflake KS, we can discuss the symbolic description of a compatible sequence 
of initial basepoints. In £)5.11 we discussed how appending an address with a 1 or 
3 essentially amounts to moving forward into a cell. If Xq is the midpoint of the 
base of the equilateral triangle, then the next point in the sequence of compatible 
initial basepoints is the midpoint of a side of Q(KSi) such that x\ and Xq are 
collinear. Symbolically, this element lies on the side si t k with address 51. Since x® 
is a midpoint of the side 51, it follows that x\ lies on the side 513. Consequently, the 
sequence of points in the compatible sequence of initial basepoints must converge to 
a point on the snowflake given by the address 513. If x\ is the basepoint of the next 
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element of the footprint of the orbit ^o(xqi Oq), then we see that Xq lies on a side 
with tt(KSo) with address 1. Since x\ is a midpoint of this side, it follows that the 
basepoint x\ of the orbit tf\(x\, 9*{) lies on the side with address 13. Consequently, 
the sequence of basepoints {xj}°°^ is converging to a point of the snowflake with 
an address 13. 

In general, if Xq is an initial basepoint of an orbit in H(KSq), then the address 
of the side s n ,k on which x^ resides is determined in a similar manner. Similarly, 
the address of the side s n ,k' containing the basepoint of the image fn(x^ , 0°) under 
the billiard map /„ associated with Q(KS n ) is also determined in such a way; recall 
that f n is a map from the phase space (KS n xS 1 )/ ~ to itself and see £|2.1l for the 
definition of the billiard map for a general polygonal billiard Vt(B). For example, 
suppose the unit interval base I of Vl(KSq) was partitioned into intervals of length 
1/3™, for some n > 0. Then suppose Xq was in one of these intervals. It follows that 
Xq would be in a similar interval of length 1/3™ if the initial direction were 7r/3. 
If this interval constitutes a ghost of some side of the prefractal approximation 
fl(KS n ), then such a ghost is removed and x\ is on some side with an address 
ending in either 1 or 3. 

Now that we have symbolically described a compatible sequence of initial base- 
points, we want to show that such a symbolic construction is entirely determined 
by the ternary representation of Xq S /. Each element of the unit interval I has a 
ternary expansion. We represent such expansions not in terms of 0, 1 and 2 but in 
terms of I, c and r, for left, center and right, respectively^ Such a coding system 
will prevent us from confusing addresses of points in / with addresses of points 
of the snowflake, but we will see that there is a very intimate connection between 
these points and a certain subset of points of the snowflake. If a ternary expansion 
of an element of the unit interval / contains an /, c or r in the (n + l)th position, 
then it follows that the element resides in either the left, middle or right third of an 
interval of length 1/3" somewhere along the base of the equilateral triangle KSq, 
respectively. 

Consider the equilateral triangle with an inscribed similar triangle sharing the 
right vertex of KSq. The side of the scaled, inscribed copy of KSq not common 
to the base and right edge of KSq constitutes a segment of a billiard ball path in 
Q(KSq), if the initial direction of the ball was ir/3. Consequently, such a leg of 
a triangle identifies a point Xq with the same ternary representation as Xq on So,fc 
such that Xq is the basepoint of fa(xQ, 7r/3), where /o is the billiard map associated 
with Q,(KSq), namely, a map from (KSq x S 1 )/ ~ to itself. For example, if xq is 
the point 3/4 in the base of Sl(KSo), then x' will also be the point 3/4, measured 
relative to the top vertex of KSq. 

We have already seen (in £14. ip that the point c (or 1/2) is the first element in 
the compatible sequence of initial basepoints giving rise to the compatible sequence 
of primary piecewise Fagnano orbits (and, therefore, each orbit has an initial angle 
of 7r/3). If Xq is some other point with a ternary expansion not equal to c, then how 
does one determine the address of the basepoint x\ of the image /„(x°,0°) from 
the ternary representation of Xq 6 II As noted above, shifting along the base of 



In §4, the use of a sequence {Pi'jli^i, with N < oo (and ji S {1,2,3} for each i > 1), 
consisting of either 0's. l's and/or 2's in order to express the decimal value of elements in the unit 
interval, was meant to allude to the ternary expansion of a number and the representation of such 
a number via an address, given here in terms of I, c and r, instead of 0's, l's and 2's. 
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the equilateral triangle results in an identical shift along any side, so long as we are 
considering a billiard ball moving in the direction 7r/3. Let us suppose that Xq € I 
has a ternary representation with c in the (n + l)th position. Then Xq will also 
have the same ternary representation and a c in the (n+ l)th position. This means 
that both points live on ghosts go,k and 5o,fc'i respectively, of Q(KSq) and will not 
be basepoints of orbits in subsequent approximations. This also implies that the 
subsequent basepoints will have addresses (in terms of the alphabet {0, 1,2,3, 4, 5}) 
such that a 1 or a 3 occurs in the {n + l)th position. Likewise, an / or an r in the 
(n + l)th position of the ternary representation for Xq dictates whether or not the 
character in the (n+ l)th position of a subsequent initial basepoint is going to have 
a 0, 2 or 4. Thus, we determine the address of sides on which initial basepoints 
reside by examining the address of Xq in the unit interval. 

5.2.1. Straightening addresses of the Koch snowflake. Let x be an elusive limit 
point of the Koch snowflake. Then x has an address that is given by a sequence 
of l's and 3's, with 0's and 2's and/or 2's and 4's interspersed, as described above. 
Without loss of generality, we may assume the address begins with a 1. At some 
point later in the address, there must be a character that is neither a 2 nor a 0. 
This character will be either a 3 or a 1. If it is a 3, we continue on to the next 
character of the address. If it is a 1, then we may switch this 1 for a 3 and until 
the next 1 or 3, switch every 4 for a while keeping the 2's unchanged. Continuing 
this process, making sure that every subsequent 1 is followed eventually by a 3 and 
not a 1, we end up with what we call the straightening of x, which we denote by 
s(x). Since s(x) has an address consisting of l's and 3's interspersed with 0's and 
2's and/or 2's and 4's such that the l's and 3's alternate, we see that this point is 
collinear with a point of / along the base of the equilateral triangle and such a line 
connecting them contains no other points of the snowflake. 

Example 5.3. Let n and k be positive integers. Consider the address of a side s U} k 
of Q(KS n ) given by a = 13123232113133100324. The straightening of this address 
is 13123212313131344120. Next, suppose the address of a particular point is given 
bya = 13123232113133100324. Then s(x) = 13123212313131344120. 

Recall from jj4.ll that M ( c (o) is dense in / and that / is the unique fixed point 
attractor of the iterated function system {fa, fa, fa} described in Equation (14.11) 
Consequently, M (the closure of M(^) in M.) is the unit interval I. Define T 
as the collection of all points of the unit interval / with finite ternary expansions 
(or equivalently, expansions terminating solely in ^'s or r's). Thus, M(^) \ T = 
I \ rJ3 Then M( c (p) \ T consists of points which have infinite ternary expansions 
not terminating solely in Z's or r's. 

If Xq were an element of T and fo{xQ, ""/•O = (^oi $o)j th en ^0 would also be a 
point with a finite ternary expansion (or an infinite ternary expansion terminating 
solely in either 0's or 2's). From this we deduce that the line segment connecting 
Xq and Xq in the direction 7r/3 would form a saddle connection in some prefractal 
billiard fl(KS n ), for some n > 0. So, when we consider allowable points from which 



28 Recall that M{V) is defined in Equation fO). 

2 ^We can then write T c (the complement of T in /) when discussing the elements of M(V)\T. 
However, so as to place more significance on the fact that such elements are limit points of M( c <?), 
we will continue using the more verbose notation. 
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to start in the direction 7r/3, we consider M( c tf) \ T. If xq is an element of this set, 
then the compatible sequence of periodic orbits {&i(x®, 0®)}°1 Q is well defined. 

Remark 5.4. Many of the elusive limit points of the Koch snowflake may be 
reached via a piecewise linear logarithmic spiral. As the name suggests, such spirals 
may be straightened to straight line segments that begin in the interior of the 
snowflake and end on the boundary of the snowflake, thus identifying an initial 
basepoint x® of an orbit of the billiard il(KSo) and the basepoint x\ of fo(x^, tt/3). 

It is clear that the addressing system and the directional exploit the inherent 
symmetry of the Koch snowflake and its prefractal approximations. In ij4.4[ we were 
implicitly using the underlying symmetry of the prefractal approximation tt(KS n ) 
as an aid in proving Theorems 14.271 and 14.291 In particular, we deferred the proof 
of Proposition 14.261 until now, because we wanted to be able to straighten a finite 
address and to explain why every basepoint of an orbit must have an identical 
ternary representation, relative to their respective sides W\ 

Proof of Proposition 14.261 Let n > and k < 3 • 4™. The point x\ on 
a side s n ^k has a finite address (where this address, in actuality, is identifying the 
side s ny k)- We know that this address can be straightened so as to identify another 
side s„.fc' of fl(KS n ) such that this side is visible from the interior of tt(KSo). In 
addition to this, one may uniquely identify x\ in terms of the address given in terms 
of {0,1,2,3,4,5} and its ternary representation given in terms of {l,c,r}. Then, 
utilizing the local and global symmetry, one can straighten x\ to an element of s n ^i . 
Without loss of generality, we may then suppose that the point on s n ^i is x° . Since 
the side s n ^< is visible from the interior of Q(KSo), so is the midpoint of s n ^k' ■ Given 
the nature of the construction of the Koch snowflake, there exist < m < n and 
h < 3 • 4 m such that s„.fc< is a subset of a side s mt h of Q(KS m ). Contained in s m .h 
is a scaled (or scaled and rotated) copy of M('if), which we denote by Mi^) m ^h- 
Since every point of M( < ^ ? ) m , ! / l ' is visible from the interior of fl(KSo), it follows that 
the midpoint of the side s m ^h is visible from the interior of U(KSq), as well. 

We can now construct a piecewise Fagnano orbit ^(j/^, 7r/3) = p&m , where the 
ternary representation of x\ and have the same ternary representation up to the 
first m many characters. This construction is precisely the content of Proposition 

EM □ 

5.3. Footprints of piecewise Fagnano orbits of J~2(KS). From $5.2[ we 
now understand how to identify points of the boundary K S and to determine the 
type of compatible sequence of periodic orbits from the initial condition (a;Q,7r/3) 
of & (x%,ir/3) in n(KS ). 

Let X = {0,1,2,3,4,5} and X n := T]-Li x = xn , the space of all words of 
length n > 1 with characters (or symbols) in the alphabet X; by default, we let 
X° := 0. Define Xoo :— YliLi x — ^ to be the space of infinite words expressed 
in terms of the elements of X. Let x 6 Xoo and n £ N. We define r„ : X^ — > X n 
by r n (x) :— x\ n , where x\ n £ X n is the finite word of length n consisting of the first 
n characters of the infinite word x. (Hence, x\ n £ X n is the truncation at level n 
of x £ Xoo.) 



"^Recall that Proposition l4.26l (from i|4.4H was used in the proof of Theorem l4.27l to calculate 
the periods of certain periodic orbits of Q(KS n )- 
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Restricting r n to KS results in a map that is well defined on the points of KS 
which are not corners of any finite approximation KSi, for any i > 0. For example, 
12 and 32 identify the same point in KS, but the words of length n given by 12. ..2 
and 32... 2 do not identify the same segments of KS n . 

Since the Koch snowflake curve KS may be viewed as the inverse limit of 
its prefractal approximations KS n , when attempting to define orbits of the Koch 
snowflake, we take as the suitable notion of limit that we have alluded to throughout 
this paper and ILaNielj to be the inverse limit 

Recall that the billiard ball path is the orbit and the footprint F n (x^) of the 
orbit & n {xP n , 7r/3) is the intersection of the orbit & n {x® n , 7r/3) with the boundary 
KS n . In H5.51 we will provide a plausibility argument as to why it may be possible 
to recast a compatible sequence of piecewise Fagnano orbits as an inverse limit 
sequence of orbits. For now, we show that the inverse limit of the footprints J-^^q) 
exists. Such a task will require us to define the proper transition maps so that 
the definition of inverse limit is satisfied. (See ^2.31 above for a brief discussion of 
inverse limits, and |Bo|, IHoYoj for further information.) In other words, we recast 
the sequence of footprints {^(^cDli^o corresponding to a compatible sequence 
of orbits as an inverse limit sequence of footprints of periodic orbits, where each 
basepoint is given in terms of its address; see ^5.11 (including Figures 121)1 |2"T1 and |2"5)) 
for a description of the addressing system. Once (and if) an appropriate billiard 
flow can be defined, we will seek in future work to show that the inverse limit of 
footprints is in some way an analogue of a Poincare section of the suitably defined 
billiard flow. 

Given m,n € N* with m < n, consider the map r nm : X n — > X m , where 
T nm (x n ) = x n \ m is the truncation of the word x n (of length n) by n — m characters 
(thereby producing a word of length m). Then T nm \KS n is a map that truncates the 
addresses of segments of KS n of length n + 1 to produce addresses of the segments 
of KSm of length m + 1. Let the map t = i mn be the 'identity map' acting on 
the unit circle S . The map l = L m n will serve to preserve the compatibility of the 
two initial conditions (x^, &m) an d i^n @n)- ^ n actuality, l(6^ ) = 8^, meaning that 
(x° m ,e° m ) is identified with (x° m ,0° n )E 

Consider XqEI with a ternary representation consisting of infinitely many c's 
and finitely many I's and r's. Recall from i)4.1l that this is a necessary and sufficient 
condition for {&i(x® , #°)}^ forming a compatible sequence of piecewise Fagnano 
orbits. Specifically, there exists a unique integer N > 1 such that &n(x%, k/3) — 

p^lf is a piecewise Fagnano orbit of H(KSn), yet for every in < N, 6 rn (x Q m , 7r/3) 
is not. Moreover, for every n > N, & n {x® n , 7r/3) is a piecewise Fagnano orbit. It 
follows that the first element x]^ in the footprint of the orbit {?n(x%,6%) has a 
finite address ending in either 13 or 31. Likewise, for each n > N, x^ has a finite 
address ending in either 13 or 31. Specifically, for every i such that N < i < n, 
(4)^0,2,4. 

Each x^ 1 has an address that can be determined from (or x^, if one prefers). 
By utilizing the local and global symmetry of the snowflake (see Figure [25]) . one 
may determine the footprint of ^ n (x°,7r/3) as a non-dynamically ordered set of 



A We note that the notion of inverse limit was used in a different context, but also in relation 
with fractals, in the recent works BeDeMiSt, RoSt, St (among others). 

^ 2 The map t maintains the compatibility of the initial basepoints and x^, by relying on 
the fact that angles are measured with respect to a fixed coordinate system. 
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points. One may then consider the inverse limit of footprints of the orbits p^f n n , 
where one forms the inverse limit by considering as our transition maps r mn (for 
m < n) the truncation of finite addresses, as defined above in this section. Letting 
J>j,(x°) be the footprint of the orbit ffnix^ , 7r/3), we write the inverse limit of the 
footprints as 



(5.1) 

hm7- n (4) = J (x*)£ N e f] -FiO^IWaftO = afr all N < m < n\ , 

{ i=N ) 

Due to the fact that this set lacks a dynamical ordering, we cannot say that 
it represents a section of any billiard flow. Therefore, we next examine how to 
formulate the footprint of an orbit of il(KS) as an inverse limit of footprints of 
orbits. 

Consider £„, the number of pairs of points preceding x„ n in the footprint T n {x G n ) 
of the orbit &n(x® n , 7r/3). One can easily see that 2£„+/3„ = k n , for some (3 n G {1, 2}. 
Then, for (a^ 4 ) G Urn J^(q;?), we have that 

(5.2) x n—l = T n-l,n( a; n n ) = T n— l,n( a 'n°" + ' 3 ") 

and x^l^ 1 is the basepoint of /^""^ 1 (a;°_ 1 ,7r/3). By the definition of a piecewise 

Fagnano orbit (see Definition 14.141) , we see that = ^-"i 1 i an d hence that 

fc n _i = Cn + 1- Consider the map F n _i ( „ given by 

(5.3) i^_ l!n (4«,^») = o 7- n _ lin x ,„_!,„ o /-*»(a£,0*»). 
The map -F n _i jn is well defined. In general, if m < n, then the map 



(5.4) F m ,„(a£»,0£") = /^ r + 1+1 o Tn _ lin x o 

constitutes the proper transition map needed to construct the inverse limit of foot- 
prints of a compatible sequence of piecewise Fagnano orbits ,8°)}f^ N . We 
denote the resulting inverse limit as follows: 



(5.5) ^(asg):=JmJi(a;J) = 

OO ^ 

{(*?«, G n = (4r ), for aU N< m < n \ . 

i=N ) 

We define a; G KS to be lim^ocX?, which is the limit (in the plane) of the 
compatible sequence of initial basepoints {x^}°^ N . 

As will be further discussed later on, we would like to think of J-(xq) as the 
'footprint' (or l Poincare section') of a 'piecewise Fagnano periodic orbiV of tt(KS), 
namely, the piecewise Fagnano periodic orbit with initial basepoint x :— linij-^oo %i 
in the direction of 7r/3. 
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5.3.1. Topological and geometric properties of the footprint .F(xq). Before dis- 
cussing our main results regarding the period, length and topological nature of 
the orbits, we establish a 1-1 correspondence between elements of M(^) and the 
collection of footprints {^(loJJiggM^)! as we h as a similar statement for the ge- 
odesic flow (and the closed geodesic paths corresponding to the unfolded billiard 
ball paths) on the associated flat surface. 

Theorem 5.5. // {^(x? , 6® )}£q * s a compatible sequence of piecewise Fagnano 
orbits, then the 'footprint' J-(xq), given by the inverse limit ^im ffj{x^, 9®), is a 
topological Cantor set {i.e., a totally disconnected and perfect compact set). More- 
over, the above inverse limit of the footprints is a self-similar Cantor set. 

Proof. Let {<%(a;?,0?)}£L o be a compatible sequence of piecewise Fagnano 
orbits. By Theorem 12.111 F{xq) is totally disconnected and compact, being the 
inverse limit of finite sets. What remains to be shown is that .F(xq) is a perfect set. 

Consider an element (x k L % 0£*)?£ o G F(xq). Fix n > 0. Then there exists 
id ,#)£o G such that 1) (x?,8?) = (yf ,0f ) for aU i < n and 2) x^ 1 ^ 

ifn+i ■ Concretely, this element (yf , <f?£ )°1 is determined from (x t * , Q i * )°^ by way 
of the local symmetry. As such, we can continue to construct a sequence of elements 
m J"(xg) that converges (with respect to either the Euclidean metric or a metric 
defined on the space of addresses) to (x k ',9-')°g a e F(x%). 

To see that such a topological Cantor set is a self-similar Cantor set, one simply 
recognizes the fact that for every footprint Fix®), there is a finite collection of IFS's 
{&x°,i}T=i) each IFS <& x o |4 consisting of two contraction mappings, and each giving 
rise to a unique fixed point attractor that is a self-similar set in its own right. 
The union of these fixed point attractors is then the self-similar footprint and one 
can easily reproduce a graphical representation of the footprint. We defer further 
explanation of this to ^5.5I 33 I □ 

Proposition 5.6. If {^{xo)} x o eM ^g^ is the collection of all footprints of piecewise 
Fagnano orbits of the Koch snowfiake billiard fl(KS) and XQ,y® £ M(f€), then 
Xq 7^ Uq if and only if JF(xq) and .F(?/q) have no elements in common. 

Proof. Suppose F{xq) ^ ^(Vo)- Recall that an address of a point can be 
straightened so as to identify a point x° g KS that is compatible with an element 
Xq (in the direction of ir/3 in the interior of KS) that is also an element of the 
unit interval base of KSq, which we have been denoting by I throughout the paper. 
This element may then be taken to be lim^oo x®, with {x®}°Z being a compatible 
sequence of initial basepoints (again, compatible in the direction of ir/3). That is, 
x° = lim^oo x® and y° = lim^oo?/?. Suppose x° ^ y°. Then there exists v > 
such that 1 1 x° — y° 1 1 > v (where we have taken || • || to be the Euclidean norm in 
W 2 ). Therefore, there exists n > such that 

(5-6) \\x° n -y°\\>v 

or 

(5-7) x° ± y° n 



More specifically, for more detailed information concerning this IFS, we refer to the dis- 
cussion immediately following the proof of Theorem 15. 121 in 35. 51 below. 
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Figure 29. A stabilizing periodic orbit of the Koch snowfiake 
billiard il(KS). The initial basepoint a: of the stabilizing orbit is 
Xi, which has the address a(x® — 5124. The initial basepoint x\ 
is compatible with the element Xq = 7/2 G /, which has a ternary 
representation of crl. 



Since x„ and Xq are compatible in the direction of 7r/3 and y^ and y^ are compatible 
in the direction of ir/3, it follows that Xq ^ y®. The converse holds since Xq ^ j/fj 
implies x a ^ y°. 

□ 

5.4. Stabilizing orbits (or ^-orbits) of Jl(KS). In §4.21 we described what 
we call an eventually constant compatible sequence of periodic orbits. Such a com- 
patible sequence of orbits was determined by the initial basepoint Xq of the first 
orbit ^q(xq, t/3) in the compatible sequence {{?i(x°, 7r/3)}°^ . The ternary rep- 
resentation of Xq is comprised of finitely many c's and infinitely many I's and r's. 
Let N be a positive integer such that for every n > N we have that & n {x Q n , 7r/3) = 
<yy(gjr.7r/3)i^ To be clear, the footprint of the orbit ^(a;^-,7r/3) consists of 
finitely many points with addresses in K S that contain finitely many l's and 3's. 
We then say that &n{x%, tt/3) is a stabilizing periodic orbit (or a ^-orbit) of 
Q(KS); see Figured 

An example of a ^-orbit of f^i^S 1 ) is one for which Xp = 7/12. This element 
is not in M(^) and the corresponding eventually constant compatible sequence of 
periodic orbits {^(x®, 7r/3)}°^ is constant after the second element of the compat- 
ible sequence of periodic orbits since the next element of the compatible sequence 
of initial basepoints corresponds to an element of the Cantor set on the side with 
address 51. 

In §4.41 Theorems 14 . 2 71 and 14.291 highlighted a relationship between the ternary 
representation of the initial basepoint Xq of &q{xq, tt/3) and the length and period 
of an orbit @ n (x® n , tt/3). The ternary representation of the element Xq also provides 
further information about the behavior of a compatible sequence of periodic orbits 
when such a compatible sequence is eventually constant. 

Theorem 5.7. Any element Xq G I with a ternary representation consisting of 
finitely many c 's and infinitely many I 's and r 's constitutes the first element in 

34 Throu ghout §4, we chose to write the ternary expansion of an element of M(V) in terms of 
the characters {0,1,2}. Because of the obvious confusion with the alphabet {0,1,2,3,4,5} that 
we use when addressing points of the Koch snowfiake KS, we use instead the characters I c and 
r to express the ternary representation of elements of M{C) and, in general, of the unit interval 
I, viewed as the base of KSo- 
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a compatible sequence of initial basepoints giving rise to a compatible sequence of 
periodic orbits that is eventually constant. 

Example 5.8. The element — crl is such an element. The eventually constant 
compatible sequence of periodic orbits stabilizes after the zeroth level approximation 
&o(xq,it/3), as depicted in Figure [231 

Further, one can determine the approximation at which the compatible se- 
quence of periodic orbits stabilizes. 

Theorem 5.9. Let Xq € I have a ternary representation consisting of finitely many 
c 's and an infinite number of I 's and r 's. Furthermore, let N be the number of left 
shifts required to eliminate all characters c from the ternary representation of the 
initial basepoint Xq of the orbit &o(xq,it/3). Then, for every n > N , we have that 
^ n «,7r/3) = n (x° n ,tt/3). 

Interestingly enough, if x a N is the initial compatible basepoint of a stabilizing 
periodic orbit, then we observe that the ternary representation Xj^ relative to the 
side sjy,fc containing x^ can be determined from the ternary representation of Xq. 

Proposition 5.10. Let x^ be the initial basepoint of a stabilizing orbit offl(KS). 
Then the ternary representation of x^ (relative to a segment of length 1/3 N ) is 
given by the Nth left-shift of the ternary representation of the point Xq. 

Example 5.11. If Xq = 5/12, then the ternary expansion of this number is clr. 
The ternary expansion of the initial basepoint x\ of the stabilizing orbit &i{x\, 7r/3) 
is given by llr. Relative to the side of length 1/3 on which this point can be found, 
the ternary expansion is Ir. 

5.5. Piecewise Fagnano 'orbits' of fi(KS). In £14.41 we saw that for every 
n > 0, there was an explicit formula for the period and length of a piecewise Fagnano 

x° 

orbit pJ^n of il(KSn)- Specifically, the length and period were determined by the 
ternary expansion of Xq. Moreover, in ^5 .31 we constructed the inverse limit of 
the Poincare sections of the orbits comprising a compatible sequence of piecewise 
Fagnano orbits {^(x?, 7r/3)}°^ . Based on the explicit formulas for the period and 
length (see Theorems 14.271 and I4.29[) , one sees that 



lim #0 n (x° n ,Tr/3) =oo 

n— too 

and, as we next explain, that 



(5.8) 



(5.9) lim \ff n (x° n ,n/3)\ 

71— >00 

is finite. 

Theorem 5.12. Let {£?i(x®, 7r/3)}°^ be a compatible sequence of piecewise Fag- 
nano orbits. Then 



(5.10) lim \0 n (x° n ,ir/3)\ 

n—too 

exists and is finite. 
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PROOF. Let Jz? be the length of the Fagnano orbit of the equilateral triangle 
billiard Q(KSq). Let {Gi(x\, 7r/3)}°^ be a compatible sequence of piecewise Fag- 
nano orbits. Then, for each n, Theorem 14.291 implies (with the notation introduced 
in Equation (14.121) ) that 



n ££ 

(5.11) \^(x° n ,n/3)\ = + Y J X[{xl\] #^-i(x°_ l! 7r/3)-. 

i=2 6 

Since, for each n > 1, \0 n (x„, 7r/3)| is bounded by \pp& n \ andlinin-^oo |ppJ^ n | = 
4«£f , it follows that 

lim |^ i (x°,7r/3)| = 2jSf + lim £ x [(sg),] tt/3)-^ 

i=a 

oo „ 

(5.12) = 2JSf + £ X [(»8)i] #^-i(x°_ 1 ,tt/3)- 

i=2 

We note that since {|(^(:e°, 7r/3)|}°^ is a monotonically increasing sequence 
(which follows from the fact that the sequence of partial sums representing this 
sequence is monotonically increasing), we deduce that the limit exists and (by 
Equation (|5.12p ) is finite. □ 

Such facts support the following plausibility argument, according to which some 
suitable limiting set can be considered as an orbit of the Koch snowflake billiard. 

Let {&i(x°, 7r/3)}^i be a compatible sequence of piecewise Fagnano orbits. 
Then there exists a least integer n > such that G n [pP n , 7r/3) is a piecewise Fagnano 
orbit and ^ > n _i(x°_ 1 ,7r/3) is not. Denote by 3 the period of ^ n _i (x°_j , 7r/3). 

Recall from the proof of Theorem 15.51 that we claimed that the footprint -F(xq) 
of a piecewise Fagnano orbit was a self-similar subset of KS. Concretely, .F(xg) is 
the finite union of (abutting) self-similar sets. Specifically, there are S many IFS's 
involved in the union describing .F(xq). That is, if for each i < 3, we have that 
$ x o ti is the operator associated with an IFS (which we denote by {$ a .o iii j}|_ 1 ) and 
^(xg)i is the unique fixed point attractor of $ x o A , then 

2 

(5.13) ■■= U^wC^oJO - H4)i- 

As a result, we have that 



j-(^) = u$ x o )i (^) i ) = u^)i- 

i—1 i—1 

For each i < 3, there exist a suitably defined non-expansive map pi that, 
together with $> x o^ tl and $3.0 ^2, constitutes a new IFS $5.0 ^, ^xg,^} ( m a 
non-standard sense). Such an IFS has an associated operator that can be used to 
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produce the piecewise Fagnano orbit p^f n ° and every subsequent piecewise Fag- 
nano orbit in the compatible sequence of piecewise Fagnano orbits. We denote the 
associated operator by 

(5-14) ^o, := Pi U $ x o tiil U Sa-g^a. 

Before showing exactly how such an operator may be utilized to produce a 
piecewise Fagnano orbit, we introduce a key notion. Let & be the Fagnano orbit 
of fl(KSo). Then, for each i < S, there is a contraction map acting on & in 
such a way that 

(5.15) v^'n ={J<t>i(.?)U0 n - l {x° n _ u x/3). 

i=i 

That is, is given by Equation (|5.15|) once a proper dynamical ordering is 

reestablished for ljf =1 ^-(J 5 ") U ff n -i(x°_ 1 ,n/3). Let &\ := where k is an 

integer meant to indicate that ^(d?) is a scaled copy of in the scaling ratio of 
1/3*. 

Finally, we have that 

i=l 

and, in general, for p > 0, 

P^i = U U ^n-l (4-1^/3)- 

j=l 

Again, a dynamical ordering must be reestablished for this equality to be true. 

Given the fact that p^ n \ v is a finite orbit with finite period, imposing such an 
ordering is not a concern. See Figure [30] and the caption therein for a detailed 
example illustrating the above discussion. 

Plausibility argument. Let {^(x?, 7r/3)}°^ be a compatible sequence of piece- 
wise Fagnano orbits and 

(5.16) p&*°° := lim (J ^„ + f(^) U ^-i(x°_ 1; tt/3), 

z=l 

where n is the least integer such that @ n {x® n , 7r/3) is a piecewise Fagnano orbit of 
ri(ifS , „) in the compatible sequence of piecewise Fagnano orbits {^(xf ,tt/3)}^ q . 
Since ^(xq) is an uncountable set and a Poincare section is at most countably infi- 
nite, it is plausible that there exists a suitably defined map / on a countably dense 
subset of J-(xq) that constitutes the Poincare section of an orbit of fl(KS) whose 
graphical representation is given by p^ x o . This may then allow us to reconstitute 
p^ x o as a true orbit of the Koch snowflake billiard ft(KS). 

In light of Theorem 15.51 and Theorem l5.12[ one would then expect this periodic 
orbit to be 'self-similar' and to have finite length. 
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Figure 30. Depicted in this figure is a piecewise Fagnano orbit 

x" 

pjF 4 of n(KS4). In addition, this orbit is an orbit in a compatible 
sequence of piecewise Fagnano orbits. We have drawn in a subset of 
the ghosts of £l(KS4) so as to see how it is one orbit is compatible 
with another. The least n such that ^„(x°, 7r/3) is a piecewise 
Fagnano orbit is n — 3. The compatible orbit ^2(^21 ?r/3) has 

z x° 

period 6. The iterated function system producing p^ 4 from pJ^ 3 
is shown. In Figure 1311 we enlarge that particular piece of the 
billiard table ft(KS4) and examine the effect of one of the IFS's 




Support for the plausibility argument. The Poincare section of p^f n \ p is finite 
and one can reestablish a dynamical ordering on Uj=i (■^fc)U^ , n -i(x°_ 1 , tt/3). 

However, the limiting set pj? 2 ^ defined by Equation (|5.16j) docs not naturally 
come equipped with such an ordering. Without additional knowledge, one cannot 
determine any dynamical ordering on p^ x » . 
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Figure 31. The effect of the IFS $ x o ti , i < 6 = #^ 2 (a^, 0%). 



Consider again the footprint J-(xq) of the piecewise Fagnano orbit with initial 
basepoint x° = limi_ > . 00 xf (and in the initial direction of 7r/3). Given how the tran- 
sition maps used in the construction of F(xq) are defined in terms of the billiard 
maps /„ defined on the prefractal billiard approximations fl(KS n ), it is reasonable 
to expect that one can recover a suitable map that would force a dynamical or- 
dering on J-(xq). This is, unfortunately, problematic, because in Theorem 15.51 we 
established that F{x < q) was a topological Cantor set, i.e., uncountable. Any attempt 
to recover a discrete set reminiscient of a Poincare section will necessarily yield an 
at most countably infinite set. Consequently, we expect to eventually be able to 
construct a suitable map defined on a countable dense subset of .F(xq) that then 
determines an appropriate analog of the Poincare section of a billiard flow for the 
Koch snowflake billiard fl(KS). With a natural dynamical ordering imposed on a 
countable dense subset of the footprint F^Xq), one may then be able to recover a 
continuous path that can be reasonably called the billiard orbit of tt(KS) with an 
initial condition (x , 7r/3). 

As mentioned, this countably dense subset would be determined by a suitably 
chosen map. The orbit itself would then be determined by how this map acts, that 
is, on what point of F(xq) such a map initially acts. As such, one may dictate that 
the initial point of departure on KS must be a point compatible with some element 
Xq G M { c &) in the direction of 7r/3 that then earns the name of a piecewise Fagnano 
orbit of Cl(KS). Thus, perhaps, the closure of the continuous orbit determined 
from a countable dense subset may correspond to the actual footprint of a piecewise 
Fagnano orbit. Moreover, if one can show that the closure of this set is the closure of 

the countable union of piecewise Fagnano orbits U^Lo P^n+pi then one has obtained 
not only a dynamical interpretation of the orbit, but also a natural geometric 
interpretation in terms of a collection of interated function systems {$ x o : i}f—i- 

5.6. Generalized piecewise Fagnano 'orbits' of f2(KS). Up until now, 
there are two types of periodic orbits of the snowflake. Those for which the initial 
basepoint of the first element of the compatible sequence is an element of Af(^), 
and those for which the initial basepoint is eventually collinear with some element 
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of a ternary Cantor set on some side s n ,k of some prefractal billiard Q(KS n ) not 
corresponding to a point with a finite ternary expansion. 

So a question naturally arises: what 'periodic orbits' (in the direction of tt/3) 
have not yet been described? Let x® £ / be an element with a ternary representation 
comprised of an infinite number of c's, Vs and r's (specifically, c's and r's, c's and 
Ps or c's, Z's and r's). For any n > 1, the compatible orbit ^ n (a;° , 7r/3) is not 
a piecewise Fagnano orbit. Nor is the compatible orbit ever an element of an 
eventually constant compatible sequence of periodic orbits. However, the existence 
of a well-defined compatible sequence of periodic orbits is not a sufficient condition 
for such a compatible sequence being recasted as an inverse limit sequence of orbits. 
Concretely, within the current framework, one cannot (at least, without further 
modifications) recast such a compatible sequence as an inverse limit sequence, let 
alone provide a plausibility argument for the existence of such an orbit, because of 
the nature of the initial basepoint Xq £ I of the orbit &o(xq,tt/3). But this does 
not preclude us from coming to terms with what exactly the limit of a compatible 
sequence of such orbits would look like. 

Conjecture 5.13 (Footprint of a generalized piecewise Fagnano orbit). Let Xq £ /. 

If 7t/3)}|2.q is a compatible sequence of generalized piecewise Fagnano orbits, 

then there exists a suitable family of transitions maps that one can use in order to 
properly formulate an inverse limit from the Poincare sections of the generalized 
piecewise Fagnano orbits. 

Conjecture 5.14 (A generalized piecewise Fagnano orbit). If one can properly 
determine the footprint of a generalized piecewise Fagnano orbit, as stated in Con- 
iecture \5.13\ then, in analogy with the support for the plausibility argument given in 
Q5.51 one should be able to recover an associated periodic orbit of the Koch snowflake 
fractal billiard fl(KS). 

6. Open Problems and Conjectures 

We close this paper by stating various additional questions, open problems 
and conjectures pointing to future research in this new field of 'fractal billiards'. 
We also offer some evidence (either theoretical or experimental) towards several 
of these conjectures. Throughout much of this discussion, we will use the present 
important and prototypical example of the Koch snowflake billiard, but naturally, 
similar questions could be asked about more general fractal billiards. 

6.1. fi(KS) eventually realized as a billiard (i.e., there may be a well- 
defined phase space) . As mentioned in the introduction, the Koch snowflake KS 
is a closed, non-rectifiable and everywhere nondifferentiable curve. This means that 
there is a priori no well-defined phase space. At the heart of the description of the 
billiard flow is the fact that we can always obtain an inward pointing vector at 
the point of collision in the boundary KS n (or any other rational billiard, for that 
matter) , modulo the conical singularities of the table. We would like to eventually 
show that a suitable analog of such a notion makes sense for the proposed Koch 
snowflake billiard fl(KS). Certainly, the absence of an inward pointing vector at 
a point x of KS would preclude us from describing an initial condition (x°, 9°) of 
an orbit. Ideally, there is an infinite collection 8 of vectors 6 of inward pointing 
directions such that at every point x £ KS and for each 9 £ 0, there exists a 
compatible sequence of initial conditions (a;°,#°)|2. such that (x, 6) £ J-(xq). 
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We therefore ask whether it is possible to find at every point of KS a collection 
of inward pointing vectors. Of course, one may then simplify the question to be 
one regarding the existence of an inward pointing vector. 

6.2. Fractal flat surfaces. In fJ31 we have shown how to construct the flat 
surfaces S(KS n ) of the associated rational billiards Q(KS n ). The genus of S(KS n ) 
is g n = 3-4" — 2; sec [MasTa, Lemma 1.2, p. 1022] for a specific formula for calculat- 
ing the genus of surfaces constructed from rational billiards. This implies that the 
sequence {g n }^Lo diverges to +00. We conjecture that there exists a surface S(KS) 
that can be appropriately viewed as a suitable limit of the flat surfaces S(KS n ) 
and such that the conjectured billiard flow on Q(KS) is dynamically equivalent to 
the geodesic flow on S(KS). If we can indeed view S(KS) as a suitable limit of 
the prefractal flat surfaces S(KS n ), then it is natural to expect that the fractal flat 
surface S(KS) should have infinite genus. 

We further conjecture that the self-similarity of the Koch snowflake curve is 
somehow reflected in the geometry of S(KS). We have gathered some experimental 
evidence towards this conjecture (see Figures 9 and 11 in |LaNiel| ). but are not 
yet ready to propose a precise formulation. 

In summary, we conjecture that S(KS) is a well-defined flat surface of infinite 
genus, that it has a fractal, self-similar structure, and that the geodesic flow on 
S(KS) should help us define (and, in fact, should be dynamically equivalent to) 
the billiard flow on the Koch snowflake table ft(KS). This is clearly a difficult and 
long-term problem that should be the focus of significant research in the future. 

Remark 6.1. When considering suitable limits of various geometric objects (such 
as flat surfaces associated with prefractal approximations to the given fractal bil- 
liards), the notion of Hausdorff-Gromov limit |Gr| of metric spaces may also be 
useful, in addition to that of inverse limit. 

Remark 6.2. We note that it would be interesting to investigate the possible con- 
nections between the types of 'fractal flat surfaces' with infinite genus conjectured 
(here and in |LaNiel|, Conjecture 4.7]) to be associated with fractal billiards, and 
the (noncommutative, adelic) 'fractal membranes' introduced by the first author in 
|La3| and viewed as (noncommutative) Riemann surfaces with infinite genus. 

6.3. Unique ergodicity of the flow. We saw in fj3]that A(KS n ) could be 
tiled by equilateral triangles A, i+ i. In the field of mathematical billiards, there is 
an important property of a particular family of fiat surfaces, called the Veech di- 
chotomy. (See [Vel-3], along with, e.g., [GuStVo,Gul,GuJul-2,HuSc,KaHa2, 
MasTa,Sm,Vo,Zo].) 

Statement 6.3 (The Veech dichotomy). For each direction 9, the geodesic flow in 
the direction 9 is either closed or uniquely ergodic. 

Since the geodesic flow and the billiard flow on a rational billiard are equivalent, 
one may rephrase the Veech dichotomy in terms of the billiard flow on a rational 
billiard Q,{B) as follows. 

Statement 6.4. For each inward pointing direction 9, the billiard flow in the di- 
rection 9 is either closed or uniquely ergodic. Specifically, independent of the choice 
of the initial basepoint x° , the path traced out by the orbit &(x°,9 ) of a billiard 
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£l(B) is either closed or uniquely ergodic in the billiard table tt(B) (as opposed to 
the billiard flow being uniquely ergodic in the phase space (B x S 1 /) 

Heuristically, a uniquely ergodic orbit of a billiard table B is one for which the 
orbit fills the tabled Likewise, a uniquely ergodic flow fills the corresponding flat 
surface. 

Theorem 6.5 (Veech's Theorem [VeT] IVe2~l IVe3) ). If the stabilizer of the flat 
surface S(B) is a uniform lattice in SL(2,M), then the Veech dichotomy holds for 
the flat surface. 

It follows from a result of Gutkin and Judge in [GuJul 2], that, for every n > 
1, the Veech group T(KS n ) (the stabilizer of the flat surface S{KS n ) in SX(2,R)) 
is commensurate with the Veech group T(KSq). Since r(A"5o) is commensurate 
with SX(2,Z), it follows from [GuJu] that, for every n > 0, the Veech dichotomy 
holds for n(KS n ). 

If a suitable surface S(KS) can be defined so that the conjectured billiard flow 
on H,(KS) is dynamically equivalent to the geodesic flow on S(KS), we conjecture 
that the Veech dichotomy will hold for the Koch snowflake billiard table ft(KS). 
To such end, we have the following theorem and definition. 

Theorem 6.6. // the flow in a direction 9 is uniquely ergodic in S(KSq), then, for 
every n > 1, an orbit in the same direction and starting at a point in S(KS n ) that 
is collinear in the direction of 9 intersects with (and stops at) at most one conic 
singularity of the surface S(KS n ), in both the past and future, thus remaining a 
uniquely ergodic direction in the surface S(KS n ). 

Proof. This easily follows from the discussion in §3 and §4. The precise 
justification, however, is somewhat lengthy, and hence the details will be provided 
in }LaNie2| . □ 

As a result, we propose the following definition. 

Definition 6.7 (Compatible sequence of uniquely ergodic orbits). If a compatible 
sequence of orbits is entirely comprised of uniquely ergodic orbits, then we say that 
it is a compatible sequence of uniquely ergodic orbits. 

Once a suitable notion of fractal flat surface has been defined and the geo- 
desic flow on such a surface has been shown to be dynamically equivalent to the 
conjectured billiard flow, a natural question to ask is whether or not the inverse 
limit (with suitably defined transition maps) of a compatible sequence of uniquely 
ergodic orbits is in fact uniquely ergodic in the billiard table Q(KS). Some work 
in progress towards this goal is being carried out in ILaNie2] . 

In summary, in the long term, one of the over-arching goals of this billiards 
project is to establish some analogue of the Veech dichotomy for the proposed 
fractal flat surface S(KS) and the associated geodesic flow. In general, we will 
attempt to establish some sufficient condition for the fractal analogue of the Veech 
dichotomy to hold for other fractal flat surfaces. 

Another important and long-term goal for the continuation of this project is 
the study of the connections between the length spectrum (i.e., the lengths of the 



■^We discuss the phase space of the billiard dynamics in i]2.1l 
It is therefore not periodic. 
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periodic orbits of the given fractal billiard) and the spectrum of the Dirichlet or 
Neumann Laplacian (or of another Hamiltonian) on the associated fractal drum 
(or 'drum with fractal boundary', as in, e.g., |Lalj . |La2| . [LaPaj . |LaN RG . 
[La-vFl §12.5], and the relevant references therein). We refer the interested reader 
to [LaNiell §4] (particularly, Open Problem 4.8) for a brief discussion of the exis- 
tence of possible Gutzwiller-type formulas [Gzl,2] (as well as [Ch,Col-2,DuGn]) 
in this context. 

In conclusion, it is hoped that the present work on the Koch snowflake billiard 
and its prefractal, rational billiard approximations, along wih its forthcoming se- 
quels (including LaNie2 ), will eventually help lay the foundations for a general 
theory of fractal billiards and their associated geodesic flows on fractal flat surfaces. 
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